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ABSTRACT 

Rich young stellar clusters produce H n regions whose expansion into the nearby molec- 
ular cloud is thought to trigger the formation of new stars. However, the importance 
of this mode of star formation is uncertain. This investigation seeks to quantify trig- 
gered star formation (TSF) in IC 1396A (a.k.a., the Elephant Trunk Nebula), a bright 
rimmed cloud (BRC) on the periphery of the nearby giant Hn region IC 1396 pro- 
duced by the Trumpler 37 cluster. X-ray selection of young stars from Chandra X-ray 
Observatory data is combined with existing optical and infrared surveys to give a 
more complete census of the TSF population. Over 250 young stars in and around IC 
1396A arc identified; this doubles the previously known population. A spatio-temporal 
gradient of stars from the IC 1396A cloud toward the primary ionizing star HD 206267 
is found. We argue that the TSF mechanism in IC 1396A is the radiation-driven im- 
plosion process persisting over several million years. Analysis of the X-ray luminosity 
and initial mass functions indicates that > 140 stars down to 0.1 M Q were formed 
by TSF. Considering other BRCs in the IC 1396 Hn region, we estimate the TSF 
contribution for the entire Hn region exceeds 14 — 25% today, and may be higher over 
the lifetime of the H II region. Such triggering on the periphery of H n regions may be 
a significant mode of star formation in the Galaxy. 



Key words: ISM: individual (Elephant Trunk Nebula, IC 1396A) - ISM: clouds - 
open clusters and associations: individual (IC 1396, Trumpler 37) - planetary systems: 
protoplanetary disks - stars: formation - stars: pre-main sequence - X-rays: stars 



1 INTRODUCTION 
1.1 Motivation 

A long-standing issue in Galactic star formation is: How 
much star formation relies on spontaneous gravitational col- 
lapse of cold clouds? vs. How much star formation relies on 
collapse "triggered" by the imposition of external forces? 
Triggered star formation (TSF) may play roles on large- 
scales from compression by Galactic spiral arms, on meso- 
scales from compression by supernova remnant (SNR) su- 
perbubbles in starburst complexes, an d on small-scales from 
comp ression by expanding H II regions (lElmegreen fc Palousl 
120071 ). It has long been recognized that star formation in 
molecular clouds can be triggered by ionization or shock 
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front s produced by nearby massive stars (|Elmegreen fc Ladal 
1 19771 ). Two major mechanisms for triggering by expand- 
ing Hn regions have been widely discussed: the Radiation 
Driven Implosion (RDI) model and the Collect and Collapse 
(C&C) model. 

In the RDI model, the expanding OB ionization 
front ablates the surface of surrounding pre-existing 
cloudlets producing bright-rimmed clouds (BRCs), often 
with cometary structures, and driving inward a compres- 
sional shock that induces star formation. Discussed since 
the 1980s, RDI now has well-developed hydrodynamical 
calculatio ns showing the viability of TSF on timescales ~ 
0.1 Mvr iKessel-Devnet fc Burkertll2003l ; iMiao et al.ll2009l : 
iBisbas et al.ll20ill . and references therein). 

Observationally, two types of bright-rimmed clouds 
(BRCs) can be distinguished: relatively massive cloudlets 
with widths of w > 0.5 — 0.7 pc, which are often found 
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to harbor more than a few stars at their tips; and thin 
(w < 0.5 pc) elonga ted pillar-like structu res with a few or no 
stars at their tips ()Chauhan et alJl201.ll ). According to the 
RDI model, this distinction may reflect different evolution- 
ary stages of BRCs with narrow pillars interpreted as the last 
vestiges of the photoevaporated and dispersed cloudlet that 
has insufficient amo unt of remaining gas to form new stars 
( Bisb as~et al.1 [201ll 'l . Alternative explanations for the thin 
elongated pill ars involve hydrodynamic instabilities in ion - 
ization fronts l|Mizuta et alj|2006l ; Iwhalen fc Normanll2008l ) 
or ionizat ion of large turbulent clo uds without pre-existing 
cloudlets l|Gritschneder et alj|201(j . and references therein). 

In most cases, molecular, infrared (IR), and Ha sur- 
veys of BRCs trace only the most recently formed stars 



^e.g.JSugitani et aljl995t lOgura et all2002i ; lUrquhart etafl 
I2009T ). While longer wavelength surveys select principally the 
youngest stars undergoing accretion from dusty disks, X-ray 
surveys efficiently select young stars due to enhanced mag- 
netic reconnecti on flaring that extends to the main sequence 
l|Feigelsonll2oTol) . In a few cases, Chandra observations have 
added the diskfree p re-main sequence (PMS) populations 
to TSF clouds (e.g.. iGetman et al.l [20071 . 120091 ). Triggered 
BRCs often show embedded mid-IR sources denoting pro- 
tostars, while Ha, JHK, and X-ray surveys reveal small 
clusters of disky and diskfree PMS stars within and in front 
of the bright rirrQ. In a few cases, spatial-age gradients in the 
stellar population are clearly seen where the youngest stars 
are embedded and older stars are ali g ned toward the ion- 
izing sources (e.g .. Qgura et al.l|2007l; IGetman et al.1 120071 ; 
llkeda et al.1 120081 ; IGetman et al.1 120*091 ). This directly sup- 
ports the RDI mechanism and implies that the existing 
clouds have been actively forming stars for more than 1 Myr. 

In the C&C model, ionization and wind shock fronts 
sweep up neutral material into a dense shell that becomes 
gravitationally unstable and fragments. Individual stars or 
entire massive clusters may form inside these fragments. 
The theory is moderately deve l oped flW hitwort h et al.1 
1 1994 iHosokawa fe Inutsukal l200fj ; I Dale et al.1 l2007al ) and 
may explain a number of spectacular examples of em- 
bedded clusters found around H n reg ions using long- 
wave l ength (mm to near-IR) obs ervations (|Deharveng et akl 
120031 ; IZavagno et ai]|2006t l2007h . They find a central ioniz- 
ing star with a roughly spherical H n region surrounded by a 
shell of dense molecular gas with dense fragments harboring 
young stellar objec ts. However, recent SPH simulations by 
IWalch et ail l|201ll) show that if an H n region expands into 
a clumpy molecular cloud then the production of a molecu- 
lar shell with dense massive fragments does not require the 
C&C mechanism. Instead, triggered star formation would 
take place through simultaneous enhancement of density and 
global radiation-driven implosion of pre-existing molecular 
clumps. 

Star formation triggered by Hn regions may be 
very common. The Spiteer/GLIMPSE survey finds 

~10 mid-IR bubbles/d eg 2 thr oughout the G alacti c 

Disk (IChurchwell et atl |200d : ISimpson et all I2012D . 
iDeharveng et al.1 J2010) show that most of these bub- 



bles enclose Hn regions ionized by 0-B2 stars, many are 
surrounded by cold dust shells and many have bright- 
rimmed dust condensations protruding inside the Hll 
region. 26% of the dust shell condensations have ultracom- 
pact Hll (UCHlls) regions and/or methanol masers, likely 
indicating triggering of massive star formation. Based on 
the cross-correlation of the positions of the mid-IR shells 
with the positions of the very bright and red young stellar 
objects (YSOs) f r om th e Red MSX Source (RMS) catalog, 
iThompson et al.1 |201ll ) estimate that the formation of 
14 — 30% of massive stars in the Milky Way could be 
triggered by expanding H 11 regions. 

Due to violent energy feedback from massive stars in 
large star- forming complexes and individual H 11 regions, the 
geometry, conditions, and history of star formation are often 
difficult to ascertain. The principal signatures of small-scale 
TSF — a star cluster in and around the cloudlet on the pe- 
riphery of an expanding HII region — may be erased through 
removal of the cloudlet by stellar winds and UV radiation, 
or the kinematic dispersal of the unbound cluster within less 
than a few million years. Thus, it is often difficult to identify 
sites of triggered star formation and to quantify the impact 
of triggering processe s. IC 1396 is a nearby (D ~ 870 pc; 
IContreras etail |2002| ) . large (~ 12 pc in radius) shell-like 
Hll region, wher e traces of recent triggered star formation 
are still evident |Kun et al.1 [2008. and references therein). 
The goal of the current study is to identify, understand, and 
quantify the triggered stellar population in the central part 
of the region associated with the largest BRC, IC 1396A 
(The Elephant Trunk Nebula), and to estimate the contri- 
bution of triggered star formation to the total stellar pop- 
ulation in both the central part and, with some additional 
assumptions, in the entire Hll region. 



1.2 The target: IC 1396A globule 

The IC 1396 Hll region, Sh 2-131 at (l,b) = (99.3°, 3.7°), 
is excited mainly by the 06.5f star HD 206267 located at 
the center of the region i n the 4 Myr old cluster TV 37 
(|Sicilia-Aguilar et al.1 120051 ) [hereafter SA05]. IC 1396 has 
a rich population of > 20 bright-rimmed and cometary 
globules seen i n silhouette against the e mission nebula 
(Sugitani et al.1 Il99ll ; IFroebrich et~ai1 120051 . and Figure [TJ 
here). Many of the clouds reside on the large molecu- 
lar shell surrounding the H 11 region, which expands at 
a speed of 5 k m s -1 with an i nferred expansion time 
around 2.5 Myr (|Patel et al.lll995l ). Sites of possible star- 
formation have b een identified in/around at least sev- 
eral globules (e . g. ISchwartz et al.1 Il99ll ; lOgura et all |2002| ; 
IFroebrich et~atl I2OO5I ) in cluding sites of substantial star- 
form a tion in IC 1396N (|Nisini et all l200ll ; IGetman et all 
120071: iBeltran et all 120091 ; IChoudhurv et alj|20ld), SFO 37 
Ikeda et al.1 \200^ U 
2 



2008 

Sicilia-Aguilar et al 



and I C 13 9 6A (iReach et al.1 |2004|; 



2009; Barcntscn et al 



20051. 120061 ; iMorales-Calderon et al.1 



201 ll ). 



1 Henceforth, the term "in front of" refers to the region between 
a globule and an ionizing star(s) in 2-dimensional projection, and 
not the region between the observer and the globule. 



Lying ~ 15' (~ 3.7 pc projected distance) west of 
HD 206267, IC 1396A is the bright-rimmed cloud closest 
on the sky to HD 206267. Its optical rim is the brightest 
of all rims and the cloud is thus likely to be the physi- 
cally closest to HP 2 06267 |Weikard et al.lll996l ). Following 
IWeikard et al.1 i| 19961 ) we designate the two brightest fea- 
tures of the rim as Rim A located in front of the > 1 pc 
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Figure X. (a) A 3° X 3° image of the emission nebula IC 1396 from the Digitized Sky Survey (DSS). The primary ionizing source of 
the region, the 06 star HD 206267, and the bright-rimmed globule IC 1396 are marked in red. The blue polygon outlines the ACIS field 
centered on the IC 1396A. (b) Close-up, combined X-ray and IR image of IC 1396A. The adaptively smoothed Chandra- ACIS image in 
the 0.5 — 8.0 keV band (blue) is superposed on the Spitzer-TRAG composite image in the 3.6 jim (green) and 8.0 fim (red) bands. The 
ACIS field is outlined in blue; the head and the trunk of the cloud are marked in red. The positions of the two brightest features of the 
ionized optical rim in front of the cloud are marked in red as Rims A and Aa. North is up, east is to the left 



diameter head of the globule and Rim Aa located further 
west at the shoulder (trunk) of the cloud (Figure [TJ. The 
part of t he cloud behind Rim Aa has its own designation: 
SFO 36 (|Sugitani et al.lll99ll ). A 0.3 pc diameter cavity in 
the molecular cloud is present in the head of the globule 
prod uced by the Herbig Ae star LkHa 349 (|Hessman et al.l 
1 19951 ). Outside this "hole", the globule is opt ically thick 
(Ay ~ 10 mag) with mole cular mass ~ 200 M jPatel et all 
ll995l ; IWeikard et al.|[l996l ). 

Early Spitzer images revealed both bright diffuse mid- 
infrared (MIR) em ission and 17 Cla ss II/I stars in the 
IC 1396A globule j Reach et all |2004 their Figure 2). A 
spectacular multi-band movie of the cloud was one of 
the first press relea s es of the Spitzer missior0. SA05 and 
ISicilia-Aguilar et all (|2006h combined new deeper and wider- 
area Spitzer photometry with optical photometry and spec- 
troscopy to reveal a rich lightly obscured population of > 200 
Ha emission and/or Li 6707 A absorption stars. The popula- 
tion is a mixture of ~ 4 Myr old members of the central Tr 37 
cluster and ~ 1 Myr stars spatially concentrated in an arc 
in front of the globule (see Figure 11 in SA05). In addition, 
~ 50 heavily obs cured IR-excess-selected ob jects were found 
inside the cloud (|Sicilia-Aguilar et al.ll2006l ). Deeper Spitzer- 
IRAC maps from the time-ser ies monitoring of the globule 
|Morales-Calderon et all 12009] ) [hereafter MC09] uncovered 
an additional dozen heavily embedded Class I/II YSOs. The 
spatial distribution of sixty five disky members in/around 
the globule is shown in Figure 16 of MC09. Recently, us- 



ing the d ata from the INT Pho tometric H-Alpha Survey 
(IPHAS). lBarentsen et~aH (|201ll ) [hereafter Bll] identified 
about two dozen low-mass high-accretion T-Tauri systems 
spatially concentrated in an arc in front of the globule (see 
Figure 18a in Bll). About half of these were newly discov- 
ered stellar members of the region. Spatial clustering of op- 
tical stars with younger ages and increasing accretion rates 
away from the ionizing star HD 206267 (SA05, Bll) led 
SA05 and Bll to propose the presence of a stellar popula- 
tion in front of the IC 1396A globule whose formation has 
been triggered by HD 206267. Thus, several high-sensitivity 
optical and mid-infrared surveys have revealed a few hun- 
dred young stars in the Tr 37 region, with the youngest 
concentrated within and around IC 1396A. 

As X-ray emission from PMS stars is based on en- 
hanced solar-type magnetic reconnection events rather than 
disk or accretion processes, X-ray selection delivers rich and 
clean samples of disk less stars missed by Ha and IR selec- 
tion (Fcigclson 2010j). Using archived X-ray Chandra grat- 



ing data. iMercer et alj |2009) identifed 22 new PMS stars 



around HD 206267, more than doubling the number of pre- 
viously reported young stars in the central 10' x 8' area 
of Tr 37. We present here the identification of > 130 pre- 
viously unknown members of the Tr 37/IC 1396A stellar 
populations, using new X-ray and optical observations with 
archived Spitzer data. The Chandra and the auxiliary Spitzer 
and optical FLWO/LAICA data are described in Iden- 
tification of X-ray sources with optical-IR counterparts and 
membership classification are considered in ij3] Classifica- 
tion of diskbearing and diskless X-ray emitting YSOs is 
provided in £|3 . 3t and the discovery of a new population of 
http://www. spitzer. caltech.edu/images/1058-ssc2003-06b-Dark-GlobulcHiMiI<C : Aa9fira IR-excess members is given in Sj4] We then es- 
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Table 1. Basic X-ray Source Properties. This table is available in its entirety in the electronic edition of the journal. A portion is shown here 
for guidance regarding its form and content. Source net counts, background counts, median energy, and apparent photometric flux are given 
for the full (0.5 — 8) keV band. Column 1: X-ray source number. Column 2: IAU designation. Columns 3-4: Right ascension and declination 
for epoch J2000.0 in degrees. Column 5: 1 cr error circle around source position. Column 6: Average off-axis angle for the two merged 
X-ray observations. Column 7: Source net counts in merged apertures and their 1 cr upper errors. Column 8: Observed background counts 
in merged apertures. Column 9: Average PSF fraction for merged observations (at 1.5 keV) enclosed within source aperture. Column 
10: Smallest of the three p-values (for the full, soft, and hard energy bands) for no-source hypothesis. Column 11: Effective exposure 
time. Columns 12: Smallest p-value for the one-sample Kolmogorov-Smirnov statistic under the no-variability null hypothesis within a 
single-observation. Column 13: p-value for the one-sample Kolmogorov-Smirnov statistic under the no-variability null hypothesis over the 
combined observations. Columns 14-15: X-ray median energy and apparent photometric flux. 
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(13) 


(14) 


(15) 


1 


213529.03+573412.3 


323.871000 


57.570087 


0.5 


11.1 


130.9 ± 12.9 


10.1 


0.90 


0.0B+00 


28.0 


1.8E-01 


1.8E-01 


2.3 


2.6E-05 


2 


213532.88+572846.9 


323.887026 


57.479703 


0.3 


10.8 


271.5 ± 17.9 


11.5 


0.91 


0.0B+00 


29.1 


3.2E-01 


3.2E-01 


1.5 


4.9E-05 


3 


213537.29+572842.0 


323.905409 


57.478345 


1.2 


10.2 


16.3 ± 6.3 


10.7 


0.91 


6.7E-07 


29.1 


6.9E-01 


6.9E-01 


1.4 


2.9E-06 


4 


213538.86+573032.7 


323.911957 


57.509103 


0.9 


9.8 


27.7 ± 7.2 


9.3 


0.90 


1.3E-12 


28.9 


8.9E-01 


8.9E-01 


2.0 


5.0E-06 


5 


213539.72+573127.5 


323.915509 


57.524320 


1.1 


9.1 


13.7 ± 5.6 


6.3 


0.89 


6.2E-07 


28.8 


6.2E-01 


6.2E-01 


3.8 


2.5E-06 


6 


213541.87+573252.8 


323.924497 


57.548001 


1.2 


9.1 


9.3 ± 5.0 


5.7 


0.89 


1.6E-04 


28.9 


9.9E-01 


9.9E-01 


2.9 


1.7E-06 


7 


213542.94+573334.8 


323.928939 


57.559685 


0.9 


9.1 


19.1 ± 6.2 


6.9 


0.90 


4.9E-12 


28.9 


6.0E-01 


6.0E-01 


1.1 


3.4E-06 


8 


213543.96+572955.9 


323.933171 


57.498871 


0.8 


8.9 


21.4 ± 6.8 


10.6 


0.90 


1.4E-08 


48.4 


2.0E-01 


2.4E-01 


1.7 


2.2E-06 


9 


213545.27+573326.5 


323.938631 


57.557375 


1.0 


8.8 


11.9 ± 5.4 


6.1 


0.90 


4.2E-07 


27.7 


7.7E-01 


7.7E-01 


1.4 


2.2E-06 


10 


213550.45+573547.5 


323.960213 


57.596543 


0.4 


9.2 


110.5 ± 11.8 


4.5 


0.90 


0.0E+00 


20.7 


2.3E-01 


2.3E-01 


1.5 


2.7E-05 



Table 2. Optical and IR Photometry of X-ray Sources. This table is available in its entirety in the electronic edition of the journal. 
A portion is shown here for guidance regarding its form and content. Column 1: X-ray source number. Columns 2-5: Optical VflwO, 
Iq FLOW-, VlaicA} and I J LAIC A magnitudes. Columns 6-9: 2MASS JHK a magnitudes, and 2MASS photometry quality and confusion- 
contamination flags. Columns 10-14: 1RAC magnitudes, and a digital flag giving photometric apertures and level of source contamination 
from nearby sources and nebular IR emission: "3" and "4" — 3-pixel (2.6") and 4-pixel (3.46") apertures with contaminating flux from a 
neighboring source of no more than 5%; "2" and "2c" — 2-pixel (1.73") apertures with possible contaminating flux of < 10% and > 10%, 
respectively. 



No. 

(1) 


v FLWO 
(mag) 
(2) 


I C,FLWO 
(mag) 
(3) 


V LAICA 
(mag) 
(4) 


I J,LAICA 
(mag) 
(5) 


.7 

(mag) 
(6) 


H 

(mag) 
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K s 
(mag) 
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F 2 
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[3.6] 
(mag) 
(10) 


[4.5] 
(mag) 
(11) 


[5.8] 
(mag) 
(12) 


[8.0] 

(mag) 
(13) 


F3 
(14) 


1 
2 
3 
4 

6 
7 
8 
9 
10 


13.44 ± 0.00 
19.50 ± 0.08 


12.37 ± 0.00 
16.41 ± 0.01 


13.57 ± 0.00 
19.48 ± 0.01 
19.44 ± 0.01 
19.36 ± 0.01 


12.01 ± 0.00 
15.51 ± 0.00 
16.48 ± 0.00 
15.71 ± 0.00 


12.64 ± 0.03 

13.65 ± 0.03 

11.49 ± 0.03 
14.30 ± 0.04 
14.22 ± 0.04 
14.61 ± 0.04 
12.86 ± 0.03 


11.82 ± 0.03 
12.95 ± 0.03 

11.18 ± 0.03 
13.44 ± 0.03 
13.38 ± 0.04 
13.87 ± 0.04 
12.28 ± 0.04 


11.56 ± 0.02 
12.71 ± 0.03 

11.07 ± 0.03 
13.07 ± 0.04 
13.17 ± 0.04 
13.66 ± 0.04 
12.11 ± 0.03 


AAAOOO 
AAAOOO 

AAAOOO 
AAAOOO 
AAAOOO 
AAAOOO 
AAAccO 


16.20 ± 0.06 
10.92 ± 0.01 
12.53 ± 0.00 
12.86 ± 0.05 
13.34 ± 0.01 


15.64 ± 0.05 
10.90 ± 0.02 
12.28 ± 0.00 
12.82 ± 0.06 
13.24 ± 0.01 


10.92 ± 0.04 
12.13 ± 0.01 
12.65 ± 0.05 
13.26 ± 0.07 


10.75 ± 0.03 
11.54 ± 0.04 


2 
3 
4 
2 
4 



timate the global properties of the Tr 37/IC 1396A stellar 
populations including age distribution spatial structure 
O, X-ray luminosity functions (XLFs) and initial mass 
functions (IMFs; SJ7J. We end with a discussion of the im- 
plications for understanding and quantifying triggered star 
formation in the region (*JH]). 



2 CHANDRA-ACIS, SPITZER-IKAC, AND 
OPTICAL OBSERVATIONS 

2.1 X-ray Data 

The current project consists of two Chandra- ACIS 
(|Weisskopf et al.ll2002l : iGarmire et aIK 2003) X-ray observa- 
tions of IC 1396A, a Guaranteed Time observation (ObsID 
No. 11807 obtained on 2010 March 31; PI Garmire) and a 
Guest Observer observation (ObsID No. 10990 obtained on 
2010 June 9; PI Getman). Both observations were pointed 



at the head of the globule but had different roll ang 
For each observation, we consider here only results arising 
from the imaging array ( ACIS-I) of four abutted 1024 x 1024 
pixel front-side illuminated charge-coupled devices (CCDs) 
covering about 17' x 17' on the sky. The S2 and S3 detectors 
in the spectroscopic array (ACIS-S) were also operational, 
but as the telescope point spread function (PSF) is consider- 
ably degraded far off-axis, the S2-S3 data are omitted from 
the analysis. The overlapping field of the two observations 
achieved 59 ks of exposure with no data losses or background 
flaring due to solar activity. 

Data reduction follows procedure s simil ar to those de- 
scribed in detail by iBroos et all |2010l , l201ll ). Level 1 event 
data were calibrated and cleaned using mostly standard 
methods and tools (CIAO version 4.3; CALDB version 
4.4.0). Event positions were adjusted to better align with 

d The aim points of the observations were 21 h 36 m 51!5, 
+57°30'51'/3 (J2000) and 21 h 36 m 47!6, +57°30'43"3 (J2000) for 
ObsIDs 11807 and 10990, respectively. The roll angles were 52? 4 
and 121? 2 for ObsIDs 11807 and 10990, respectively. 
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Table 3. Derived Properties, Stellar Identifications, Membership, and Classification. This table is available in its entirety in the electronic 
edition of the journal. A portion is shown here for guidance regarding its form and content. Column 1: X-ray source number. Columns 2-3: 
X-ray net counts and median energy for the full (0.5 — 8.0) keV band. Columns 4-5: X-ray column density, intrinsic luminosity, and their 
errors (summed in quadrature statistical and systematic errors). Luminosities are derived assuming a distance of 870 pc and are given for 
the full (0.5 — 8.0) ke V band. Columns 6 -7: Age estimates derived from the optical V vs. V — Iq an d V vs. V — Ij color magnitude diagrams 
using PMS models of 1 iess et al. and assuming a constant source extinction for optical PMS candidates of Ay = 1.56 mag (SA05, 

Bll). Column 8: Apparent SED slope from IRAC photometry with 1<t error. Column 9: Number of IRAC bands, from which the SED 
slope was derived. Column 10: 1MASS identifier. Column 11: Source identifier from our analysis of the IRAC data; X-ray sources outside 
the field of view of the IRAC data analyzed here are labeled as "OutOflracFOV" . Columns 12-14: Stellar counterparts from SA05, Bll, 
and MC09. Column 15: This positional flag indicates if an X-ray source lies projected against the globule ("1") or not ("0"). Column 16: 
Indicates membership and YSO class: "EXG" — possible extragalactic contaminant; "FRG" — possible foreground contaminant or YSO; 
"DSK" and "NOD" — disky and diskless YSOs, respectively; "UNCI" and "UNC2" — objects of the uncertain class. "UNCI" are sources 
with no IR counterparts; these are most likely extragalactic objects but some can be background stars. "UNC2" have weak registered and 
unregistered IR counterparts or lie close to bright IR objects; these could be YSOs or field stars or extragalactic objects. 
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the Naval Obs ervatory Merged Astr ometric Dataset cata- 
log (NOMAD: IZacharias et al.l l2004h . The effective point- 
spread function was slightly improved by removing a pixel 
randomization that was historically added to ACIS event 
positions, and by applyin g a sub-pixel position estimation 
algorithm (|Li et alj |2004| ). A custom bad pixel table was 
used to slightly increase point source sensitivity. Instrumen- 
tal background events were identified and removed using an 
aggressive algorithm when searching for sources and using an 
algorithm with few false positives when extracting sources 
l|Broos et al.ll2010l . their Figure 1). 

Detection of point sources is performed usin g both 
the wavelet transform method (| Freeman et alj |2002| ) 
and the Lucy-Rich ardson image reconstruction algorithm 
l|Broos et al.l l201p| . their §4.2); the latter is better suited 
for resolving closely spaced sources. Updated position es- 
timates, photometry, and spectra were obtained for candi- 
date sources using the A CIS Extract (AE) software packagtQ 
jBroos et al.ll2010l . l2012l ). 

Some of the inferred properties, such as X-ray net 
counts (NC), background counts (BC), median energy 
(ME), and apparent flux (Fx) for the full (0.5 - 8) keV 
energy ba nd, as well as p-v alue for no-source hypothesis 
(PhNoSrc: iBroos et al"1l2010l , their §4.3 and Appendix B), 
are provided in Table[T] Similar t o the Chandra cata log of X- 
ray sources in the Carina Nebula |Broos et al . 2011), our list 
of candidate sources in IC 1396A is trimmed to omit sources 
with fewer than 3 total source counts (NC + BC < 3) and 
probability for being a background fluctuation greater than 
1% (P b NoSrc > 0.01). 

Our final catalog comprises 415 X-ray sources (Table 
[1]). We expect that roughly half of these sources are extra- 
galactic with extremely optically-faint counterparts ( i)3.2[) . 



4 The ACIS Extract software package and User's Guide are avail- 
able at http://www.astro.psu.edu/xray/acis/acis_analysis.html . 



and the rest are young stars associated with the Tr 37 and 
IC 1396A star forming regions. 

2.2 IRAC Data 

To obtain MIR photometry for X-ray objects and to identify 
and measure MIR photometry for additional non- Chandra 
disky stars that were missed in previous studies o f the region, 
we h ave reduced the archived Spitzer-IKAC ijFazio et al.l 
|2004| ) data from a Guaranteed Time observation of the 
program #37 (Fazio, PI). The same observation has been 
presented in SA06 and lSicilia-Aguilar et al.l (|201ll ) [hereafter 
SA11]. The SA06 data had been reduced using old calibra- 
tion and flat fielding information, and the photometry for 
sources with highly non-uniform background was addition- 
ally affected by the choice of large source extraction aper- 
tures. Given the improvement of the pipeline and calibration 
since then, SA11 have re-reduced the data and redone the 
photometry for a selected sample of YSOs. Since both SA06 
and SA11 catalogs lack photometry measurements for the 
bulk of the X-ray sources, reanalysis of the data is in order. 

The observation was obtained on 2003 December 20 
with the IRAC detector in all four IRAC channels, 3.6, 4.5, 
5.8, and 8.0 /im. Two adjacent fields subtending ~ 37' x 42' 
in channel pairs 3.6/5.8 fim and 4.5/8.0 fim were centered 
on Tr 37. These fields were imaged with an 8 x 9 mosaic of 
adjacent positions separated by 280"; individual exposures 
were taken at 5 dithered positions for each of the 72 cells 
of the mosaic. The images were taken in the high-dynamic- 
range mode with 30 and 1.2 second frame time exposures. To 
reduce unnecessary data processing we have analyzed only a 
portion of the original data that encompasses the Chandra- 
ACIS field with the coverage of ~ 19' x 19' area in all four 
channels centered on the Rim A of the IC 1396A globule. 
This covers 93% of the ACIS field omitting its north-western 
and south-western edges. 

Basic Calibrated Data (BCD) products from the Spitzer 
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Figure 2. Entire IRAC source sample. [3.6] band. The photometry from 4 pixel wide background annuli versus the photometry from 
1 pixel wide background annuli. The top 4 panels show the photometric measurements for the source samples with different source 
extraction apertures: 4-pix (orange), 3- pix (blue), 2-pix (green), and 2-pix with possible contaminating flux from a neighboring source 
of > 10% (red). The bottom panel is a histogram representation of the photometric differences for all 4 source samples with different 
source extraction apertures. The median, mean, and standard deviation of the magnitude differenceas are: 0.069, 0.094, and 0.133 mag 
for the "2c" extractions; 0.02, 0.034, and 0.099 mag for the 2-pix extractions; 0.007, 0.011, and 0.077 mag for the 3-pix extractions; and 
-0.008, -0.008, and 0.06 mag for the 4-pix extractions. 
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Science Center's IRAC pipeline version S18.18.0 were auto- 
matically treated with the WCSmosaic IDL package devel- 
oped by R. Gutermuth from the IRAC instrumental team. 
Starting with BCD data products, the package mosaics indi- 
vidual exposures while treating bright source artifacts, cos- 
mic ray rejection, distortio n correction, subpixel o ffsetting, 
and background matching l|Gutermuth et all 120081 '). We se- 
lected a plate scale of 0.86" for the reduced IRAC mosaics, 
which is the native scale divided by y/2. 

Aperture photometry of IRAC sources was obtained 
using the IRAF task PHOT. Circular apertures with ra- 
dius of 2, 3, and 4 pixels (1.72", 2.58", 3.44") are assigned, 
with larger apertures favored for uncrowded sources and 
smaller apertures favored for crowded sources. Using the 
IRAC point-spread function (PSF) for the [3.6] band resam- 
pled to a plate scale of 0.86", pairs of sources are simulated 
with wide ranges of source separations, orientations, and flux 
ratios to derive the flux contribution from a nearby source 
within 2, 3, and 4 pixel apertures as a function of separa- 
tion angle and flux ratio (Kuhn et al. in preparation). For 
the real sources, their photometry is derived using 2, 3, and 
4 pixel source apertures, and we report photometry from 
the largest aperture for which the expected contamination 
from a nearby source is less than 5%. A source aperture ra- 
dius of 4 pixels was used for most of the X-ray sources and 
non- Chandra selected sources of interest O- The cases for 
which the expected contamination is > 10% in the smallest 
2 pixel aperture are flagged in Table [2] as "2c" . 

The background in the region is highly non-uniform due 
to the contribution by the PSF cores and wings of multiple 
neighboring sources often superimposed on top of the spa- 
tially variable diffuse emission. For all IRAC sources the 
aperture photometry was performed for two different sizes 
of sky annuli, 4 pixel and 1 pixel wide sky annuli adjoin- 
ing source extraction regions. The former, closer in size to 
the "standard" sky annuli recommended in the IRAC In- 
strument Handbook, collects higher numbers of sky pixels 
offering smaller statistical uncertainties but with the dan- 
ger of increased systematic error due to sampling wrong 
background. Th e latter is chosen i n accordance w i th th e 
strategy of e.g . . lLada et alJ (|2006l ): iLuhman et all l|2008h : 
iGetman et guTl|2009h and K. Luhman(private communica- 
tion) who attempt to sample the actual sky present within 
the source extraction region. For background selection, the 
optimal balance between proximity and number of back- 
ground pixels is not obvious. 

We adopt zero point magnitudes (ZP) of 19.670, 18.921, 
16.855, and 17.394 in the 3.6, 4.5, 5.8, and 8.0 /mi bands , 
where M = -2.5 \og{DN/sec) + ZP (|Reach et al.ll2005h . 
Since our plate scale and the extraction apertures have cus- 
tom sizes, the aperture correction val ues are calculated her e 
using the current data as desribed bv lLuhman et alJ i|2010h . 
The calculations are based on the comparison of the photo- 
metric measurements for the relatively isolated, bright, and 
unsaturated sources obtained for the source and sky regions 
of different sizes including large 14 pixel (~ 10 native pixel) 
extraction regions that are sim il ar in size to those used in 
the calibrations of iReach et al.l (|2005l ). The total aperture 
corrections applied to our measurements in the case of 4-pix 
background annulus are: 0.213, 0.184, 0.218, 0.416 mag in 
the case of 4-pixel source aperture; 0.441, 0.367, 0.698, 0.949 
mag in the case of 3-pixel source aperture; and 0.829, 0.955, 



1.319, 1.105 mag in the case of 2-pixel source aperture in 
the 3.6, 4.5, 5.8, and 8.0 ^m bands, respectively. The total 
aperture corrections applied to our measurements in the case 
of 1-pix background annulus are: 0.174, 0.171, 0.161, 0.180 
mag in the case of 4-pixel aperture; 0.387, 0.304, 0.496, 0.653 
mag in the case of 3-pixel aperture; and 0.637, 0.714, 0.976, 
0.868 mag in the case of 2-pixel aperture in the 3.6, 4.5, 5.8, 
and 8.0 /im bands, respectively. The reported photometric 
errors include statistical errors in the source and background 
emission and a 2% uncertainty in the calibration of IRAC 
l|Reach et all 120051 ). For further analyses we retain IRAC 
sources with signal-to-noise > 3 in both the [3.6] and [4.5] 
bands, and the photometry is reported if the signal-to-noise 
> 3. 

For all extracted IRAC sources we compare their photo- 
metric measurements performed with 1-pix and 4-pix back- 
ground annuli. Figure [2] exemplifies such a comparison for 
the [3.6] band. For the sources with source apertures 4-pix 
(orange), 3-pix (blue), and 2-pix (green) the distributions 
of the magnitude differences typically have small biases of 

< 0.02 — 0.05 mag and dispersions of < 0.04 — 0.2 mag 
increasing towards the fainter magnitudes and longer wave- 
lengths. Both the effect of the finite number of background 
pixels and the possible physical mistakes the two methods 
are making could be the cause of the magnitude biases and 
spreads. In the case of the most crowded sources, with the 
apertures "2c" (red), the bias and dispersion are higher, 

< 0.02 — 0.08 mag and < 0.1 — 0.4 mag, respectively. The 
nature of this bias is unclear. For example, if a neighbor lies 
within the 1 pixel aperture, then the larger 4 pixel aper- 
ture would dilute the light from this neighbor, shifting the 
photometry in the observed direction. For the X-ray sub- 
sample of "2c" sources the bias appears only in the range of 
[3.6] > 16 mag characteristic of objects unrelated to this rel- 
atively close star forming region f H3.2p . Performing the disk 
classification based on spectral energy distributions (SEDs) 
as described in i]3.3l we find that the choice of the back- 
ground aperture appears to have a tiny effect on SED shapes 
and no discernible effect on disk classification for all X-ray 
YSOs. Since there is no evident preference for either back- 
ground annulus size, we publish the photometry from the 
method with the 1 pixel wide background annuli, which pro- 
duces conservatively large statistical errors. Derived IRAC 
photometric measurements for 242 X-ray sources are given 
in Table [2] 

For the objects in common with MC09 we compare 
our photometric measurements with the time-averaged mea- 
surements of MC09 that are based on new deeper Spitzer- 
IRAC exposures. Figure [3] exemplifies such a comparison for 
the [3.6] band. The source sample here is relatively bright 
([3.6] < 14 mag), so the magnitude differences are not af- 
fected by the choice of our background apertures. The mag- 
nitude differences have small biases of 0.02 — 0.03 mag and 
dispersions of 0.22 — 0.25 mag. About 60% of the sources 
considered here are MIR variable. Since the sources are rel- 
atively bright, variability rather than photometric uncer- 
tainty is likely the major cause of the magnitude spread. 
For the 28 X-ray YSOs with MC09 counterparts, our SED- 
based disk classification agrees with that of MC09 for all but 
one source (JO}. 

Spitzer-MlPS data are also available (PID 58; PI Rieke). 
The use of the MIPS data in the current work could poten- 
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Figure 3. Photometry comparison with the time-averaged measurements of MC09. [3.6]-band. Left (right) panels present the photometric 
measurements for our 1-pix (4-pix) background methods versus MC09 photometry. The top and bottom panels give a scatter-plot and 
a histogram versions, respectively. The source samples are color-coded according to our choice of different source extraction apertures: 
4-pix (orange), 3-pix (blue), and 2-pix (green). The MfR- variable sources from MC09 are additionally marked by the red squares. On 
panel (c) the median, mean, and standard deviation of the magnitude differenceas are: -0.102, -0.207, and 0.341 mag for the 2-pix 
extractions; 0.062, 0.057, and 0.203 mag for the 3-pix extractions; 0.056, 0.064, and 0.153 mag for the 4-pix extractions; and 0.026, 0.003, 
and 0.237 mag for the entire sample. On panel (d) the median, mean, and standard deviation of the magnitude differenceas are: -0.082, 
-0.176, and 0.316 mag for the 2-pix extractions; -0.058, -0.007, and 0.224 mag for the 3-pix extractions; 0.056, 0.07, and 0.141 mag for 
the 4-pix extractions; and 0.024, 0.006, and 0.222 mag for the entire sample. 



tially lead to revised disk classes for some of the sources — 
from diskless classified by shorter wavelengths to Transition 
Disks (TDs), systems with inner disk holes. However, due to 
the reduced spatial resolution and sensitivity of the MIPS 
data such a classification would be biased towards brighter, 
less crowded objects located mainly outside the globule (the 
region with the brightest diffuse emission). For simplicity, 
in this work we prefer to omit the usage of the MIPS data 
and to classify any TD candidate as diskless ( £I3.3[) . This 
should not affect the goals of the paper, identification and 
characterization of triggered stellar populations. 



2.3 FLWO and Calar Alto Data 

The optical observations of the Tr 37 an d IC 1396A and their 
data analyses are described in detail bv lSicilia-Aguilar et al.l 
(|2004L I2005L l2010h . 

The UVRcIc observations were carried on with the 
1.2 m telescope at the Fred Lawrence Whipple Observatory 
(FLWO), using the 4Shooter C CD array, between Septem- 
ber 2 000 and September 2002 (|Sicilia-Aguilar et al.ll20ol . 
2005). The 4Shooter is an array of four CCDs, covering a 



square of 25' on the side. Two 4Shooter fields were taken to 
cover ~ 45' x 25' area centered on HD 206267. The FLWO 
fields contain the whole ACIS field, except for a small gap in 
between the 4 CCDs of the 4Shooter. All but a few Chandra 
stars were observed in September 2000. 

The UVRjIj observations of Tr 37/IC 1396A were ob- 
tained in service mode during three nights in 2007 June 9-11, 
using the Wide-Field Camera L AIC A mounted on the 3.5 m 
telescope in Calar Alto, Spain (jSicilia-Aguilar et al.ll201(ih . 
LAICA is a 2 x 2 mosaic of four CCDs, each covering a 
15.3' x 15.3' field of view with a large gap of 15.3' x 15.3' in 
between. The project combines four LAICA pointings cov- 
ering ~ 45' x 45' area around HD 206267, including nearly 
the entire ACIS field. 

After applying standard procedures for bias and flat 
field cor rections, aperture photometry, a nd zero point cali- 
bration ijSicilia-Aguilar et al.|[20oi . l2010h . we obtain optical 
photometry for 181 X-ray sources. The resulted Vic and 
VI j photometric measurements are given in Table [2] 
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3 INFRARED AND OPTICAL 

COUNTERPARTS TO X-RAY SOURCES 

3.1 Identifications 

Positions of our X-ray sources were compared with source 
positions from the three optical catalogs (SA05/FLWO, Bll, 
and LAIC A), the NIR 2MASS catalog, and the two MIR 
Spitzer catalogs (MC09 and our custom catalog from H2.2p . 
SA05, Bll, LAICA, and 2MASS catalogs cover the full 
ACIS field. Our custom MIR catalog overlaps 93% of the 
ACIS field omitting its north-western and south-western 
edges. MC09 catalog covers the inner region of the ACIS 
field near/inside the globule; this amounts to ~ 20% of 
the ACIS field. MC09 catalog inclu des all diskbearing stars 
near/inside the globule rep orted bv lReach et al.l (|2004T ) and 
ISicilia-Aguilar et all l|200ej) . 

An automated cross-correlation between the Chandra 
source positions and optical-IR source positions was made 
using a search radius of 2" within ~ 6' of the ACIS field 
center, and a search radius of 3.5" in the outer regions of 
the ACIS field where X-ray source positions are more uncer- 
tain due to the deterioration of the Chandra telescope point 
spread function. This was followed by a careful visual ex- 
amination of each source in both bands to remove dubious 
sources and associations. 

Optical FLWO/LAICA (323]), NIR 2MASS, and MIR 
IRAC f £|2.2[) photometry for X-ray sources is given in Table 
Individual counterparts from 2MASS, our custom IRAC 
( fl2~2l) , SA05, Bll, and MC09 catalogs are given in Table[3] 

Chandra detects 27 out of 69 (39%) of the MC09 
infrared-excess objects and 45 out of 76 (59%) of the 
B11/SA05 Ha emitters and Li absorbers. Most of the optical 
YSOs undetected in X-rays are diskbearing stars. The X-ray 
detection efficiency of disk bearing stars is known to be lower 
than that of diskless stars (|Getman et al . 2009; Stelze r~et alj 
l201ll . and Figure [S] here) . To compensate for this effect, in 
section i|4]we introduce a catalog of additional non- Chandra 
IR-excess selected disky stars within the ACIS field and be- 
yond. This allows a recovery of 83% of Ha emitters and Li 
absorbers. Several optical stars located at the very edge of 
the ACIS field, where the ACIS sensitivity degrades, remain 
unrecovered. 



3.2 Pre-Main-Sequence Census 

X-ray surveys of star-forming regions suffer contamination 
by extragalactic sources, mainly quasars and other active 
galactic nuclei, which can be seen even through the Galac- 
tic plane as faint, absorbed X-ray sources. Additional con- 
tamination arises from foreground and background Galactic 
stars, mainly main-s equence stars and some types of giants 
l|Getman et al.ll201ll . and references therein). 

We perform detailed simulations for extragalactic and 
Galactic X-ray contamination populations expected in the 
direction of the IC 1396A ACIS field, (l,b) = (99.1°, +3.9°). 
The methodo l ogy for such simulations is detailed in 
iGetman et alj l|201lf ). The simulations take into considera- 
tion a variety of factors involvi ng a Galactic population syn- 
thesis model l|Robin et al.ll2003h . stellar (AGN) X-ray lumi- 
nosity (flux) functions, Chandra telescope response, source 
detection methodology, and possible spatial variations in 



the X-ray background and absorption through molecular 
clouds. For the IC 1396A cloud, we construct its extinc- 
tion map based on dust redenning of 2MASS sources exclud- 
ing possible foreground sources, all previously known YSOs, 
and 2MASS sources with X-ra y counterparts (for de tailed 
methodology see Appendix B in lSchneider et al"1l201ll ). The 
derived extinction varies from Av < 2 mag outside the glob- 
ule to 3 < Av < 12 mag through the head of the main glob- 
ule (behind Rim A) and along the edge of SFO 36 cloudlet 
adjacent to Rim Aa (Figure[T}. Our simulations predict that 
~ 20 — 30 foreground X-ray stars with X-ray median ener- 
gies of ME < 1 keV, ~ 20 — 30 background stars with 
1 < ME < 2 keV, and ~ 110 — 150 extragalactic objects 
with ME > 2 keV are present in the ACIS field. 

Based on presence/absence of IR counterparts, source 
spatial distributions, X-ray counts, X-ray median energies 
(surrogates for X-ray absorption), and the predicted num- 
ber of contaminants, the X-ray objects detected in our Chan- 
dra observation can be separated into the following source 
classes (see Table [3] and Figure |4)|: 

A. Bright extragalactic objects ("EXG"). 

B. Uncertain class 1 (labeled as "UNCI") mostly com- 
posed of faint extragalactic objects. 

C. Uncertain class 2 ("UNC2"), which is possibly a mix- 
ture of contaminants and some low-mass YSOs. 

D. A mixture of foreground stars and possibly some YSOs 
("FRG"). 

E. Diskless and disky young members of the region 
("NOD" and "DSK"). 

The 29 hard X-ray sources (ME > 2 keV; magenta in 
Figure H (a)) from the class "EXG" lack 2MASS but have 
dim MIR counterparts ([3.6] > 16 mag, [4.5] > 14.5 mag). 
These are most likely bright extragalactic objects, because 
they are clear outliers from the stellar locus on the [3.6] 
versus X-ray net count diagram (Figure U (b)), and their 
[3.6] and [4.5] magnitudes and X-ray fluxes are fully consis- 
tent with those of the Chandra- Spitzer AGN sample from 
the Serendipitous Extragalactic X-ray Source Identification 
(SEXSI) program (lEckart et al.ll2010h . 

The class "UNCI" is composed of 162 X-ray sources. We 
believe that most of these sources are extragalactic objects 
because they lack IR counterparts, are generally X-ray faint 
and hard (NC < 4 - 10 counts and ME > 2 keV; Figure 
2] (a)), and their number is comparable to the number of 
simulated AGNs. Their spatial distribution is also similar to 
that of the simulated extragalactic objects: the weakest X- 
ray sources are concentrated towards the central part of the 
ACIS field while brighter sources are uniformly distributed 
across the field (black points in Figure|4](c)-(d)). The greater 
density of very weak X-ray sources on-axis than off-axis is 
due to variation in detection completeness with off-axis angle 
l|Broos et al.ll201lh . The detected X-ray YSOs are much less 
prone to such an effect because only a handful of them are 
very weak sources (red and blue points in Figure |4]). 

The 28 objects of the class "UNC2" do not have cat- 
aloged IR counterparts, but our visual inspection of the 
2MASS and Spitzer-IRAC images suggests that some might 
have very dim uncatalogued counterparts and some lie near 
bright IR objects. We suggest that this class is a mixture of 
heterogeneous contaminants and low-mass YSOs. 

The remaining X-ray sources with IR counterparts are 
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Figure 4. Properties of different classes of X-ray sources. Source classes are color-coded: the uncertain class 1 (UNCI), which is largely 
comprised of extragalactic objects, is in black; the uncertain class 2 (UNC2), which is possibly a mixture of contaminants and some 
low-mass YSOs, is in cyan; the EXG class, which is composed of IR- X-ray bright extragalactic objects, is in magenta; the FRG class that 
likely contains foreground stars and some YSOs is in green; and the NOD/DSK classes are composed of young dikless/disky members 
shown in red/blue, (a) X-ray net counts versus X-ray median energy. The dashed line marks the 4 net count level. One of the most 
reddest Class I protostars found within the globule {Chandra source #32 = IC1396A-7) is labeled, (b) 3.6 /xm-band magnitude versus 
X-ray net counts, (c)-(d) Spatial distribution of different classes of X-ray sources superposed on the Chandra exposure map shown in 
inverted colors with logarithmic scales. The distributions are stratified by X-ray net counts, NC < 4 counts (c) and NC > 4 counts 
(d). The green polygon outlines the Chandra FOV. The red contour demarcates the globule, tracing 8.0 fim emission from hot dust and 
PAHs. 



further divided into 21 possible foreground stars (class 
"FRG") and 175 likely YSO members of the region (classes 
"NOD" and "DSK"). The "FRG" class is composed of stars 
with low median energies and NIR colors (green in Figure 
[4] (a) and Figure [8] (a)) that are not listed in the optical 
catalog of young stars from SA05. 

The application of the MIR source classification scheme 
of iGutermuth et alj (2009) to the X-ray YSO sample yields 
zero extragalactic candidates. The scheme can not be ap- 
plied to the "EXG", "UNCI", and "UNC2" objects, since 
many of them lack MIR counterparts and all lack [5.8] and 
[8.0] band magnitudes. 

Only the 175 likely YSO members (classes "NOD" and 
"DSK") are used in our analyses of the IC 1396A stellar 
population. In £13.31 these sources are further separated into 
124 diskless and 51 disky X-ray YSOs. The J versus J — 



H color-magnitude diagram in Figure [5] shows that the X- 
ray catalog is sensitive to diskless and disky stars down to 
J > 15 mag and J > 14 mag, respectively. X-ray YSOs 
outside the globule have typical source extinctions of Av < 
2 mag, while inside the globule stars can be subject to higher 
extinction of up to Ay = 7—10 mag. The optical catalog 
of SA05 within the ACIS field (green o) is not sensitive to 
stars with J > 14 mag or highly-absorbed stars. Figure [S] 
shows that at J > 14 mag the X-ray disky sample is well 
complemented by the MC09 catalog of disky stars near and 
inside the globule as well as our custom IRAC catalog of 
disky stars outside the globule (Q . 

In the ACIS field, previous studies had located 76 Ha 
emitters and Li absorbers (SA05, Bll) and/or 69 mid- 
infrared excess stars (MC09) representing Class I/II YSOs. 
The Chandra source list (Tables IH3J1 includes 61 of these 
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Figure 5. NIR color-magnitude diagram for all Chandra YSOs with available NIR photometry. Diskless and disky X-ray stars are 
marked as red and blue X, respectively. Sources projected against the globule are further outlined by green □. X-ray sources w ith known 
optical counterparts from SA05 are additionally marked by green o. 1 and 4 Myr PMS isochrones from [siess et all {2000) at 870 pc 
distance are shown as solid black and magenta lines, respectively. Reddening vectors of Ay = 5 mag for various masses are indicated by 
dashed lines. Some previously known intermediate and high-mass stars are labeled. 




Figure 6. NIR color-magnitude diagram for all Chandra+Spitzer disky stars with available NIR photometry identified within the ACIS 
field. X-ray stars are marked by blue X. Disky stars from MC09 are marked by green o. Additional non- Chandra and non-MC09 disky 
stars (from Q within the AC IS field ar e mark ed as cyan □. Stars projected against the globule are further outlined by green □. 1 and 
4 Myr PMS isochrones from lsTess et all I teooch at 870 pc distance are shown as solid black and magenta lines, respectively. Reddening 
vectors of Ay = 5 mag for various masses are indicated by dashed lines. 
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Figure 7. IR SEDs for 14 new X-ray disky stars that are not listed in the previous YSO catalogs. Plus, a diskfree SED for the highly 
absorbed X-ray source # 235. JHK S (diamond) and IRAC-band (square) flux point s with usually small errors. The dashed and solid 
lines give the original and (de) reddened IC 348 median SED from lLada et al.l {200g) fitted to the data. The top two lines of the panel 
legends give information on stellar mass and source extinction estimates derived from the c olor-magnitude J v ersus J — H diagram. The 
third line gives the spectral class of the corresponding IC 348 median SED template from lLada et al.l J200&-) and reddening applied to 
the original template SED to fit the observed IC 1396A source SED. 
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previously known members, and discovers 114 new mostly 
diskless member^]. 



3.3 Pre-Main-Sequence Disk Classification of 
X-ray stars 

We base the evolutionary classification of X-ray emitting 
YSOs in the 1C 1396A region on a comparison of their IR 
spectral energy distributions (SEDs) with the SEDs of YSOs 
in the well-studied IC 3 48 cluster in the Perseus molecular 
cloud l)Lada et al.ll2006T h Specifically, we compare the ob- 
served IC 1396A SEDs in 2AfA5'S+IRAC IR bands to the 
(de)reddened median SED templates of IC 348 PMS stellar 
photospheres. This is ou r primary class i ficatio n method and 
is described in detail in lGetman et al.l (120081 . their §3). We 
further confirm these classifications through calculations of 
the observed SED slope ao, and locations in IRAC color- 
color diagrams. Results of our disk classification are given 
in Table El 

In the SED-based classification, rough estimates of in- 
dividual stellar masses of the IC 1396A X-ray YSOs are 
obtained from the J versus J — H color-magnitude diagram 
(Figure [5} using a set of trial age values (1, 2, 3, 4, and 
10 Myr). For each trial age, an IC 348 template SED at 
the corresponding spectral sub-class is assigned to each of 
the X-ray sources. The IC 348 template is (de)reddened and 
normalized to match both J and //-band fluxes of the ob- 
served IC 1396 A source's SED. The reddening applied to the 
original IC 348 templates to fit the observed IC 1396A SEDs 
(AAy in Figure [TJl summed with the median source extinc- 
tion for the diskl ess YSO sample in IC 348 (Av =1.6 mag 
from Table 2 in lLada et al l teoOrJ n gives an estimate of 
the source extinction from the SED fitting, Av,sed- We 
find that Av,sed is in a reasonable agreement with the 
source extinction derived from the J versus J — H color- 
magnitude diagram, Av,cmd- For example, for the trial age 
of 2 Myr the median, mean, and standard deviation of the 
Av,sed — Av,cmd difference are —0.1, —0.1, and 0.6 mag, 
respectively. 

Visual examination of deviation of the observed IC 
1396A SED from the corresponding matched IC 348 pho- 
tospheric template allows disk classification of IC 1396A X- 
ray YSOs. While inferred values for individual stellar masses 
strongly depend on the choice of age, the final classification 
discriminating diskless and disky stars is not age sensitive. 
With regard to the 28 X-ray YSOs with MC09 counterparts, 
our disk classification agrees with that of MC09 for all but 
one (Chandra source #235 = MC09 source IC1396A-46) 
sources. Disky SEDs for 14 X-ray YSOs that do not have 
counterparts in any previous stellar catalog and a diskfree 
SED for the highly absorbed X-ray YSO #235 are shown in 
Figure [7] 

We also performed least-squares linear fits to the 
log(A/\) values in the four IRAC wavelength band to ob- 
tain the observed (not de-reddened) SED spectral index, an. 



5 Notice on Figure[4](a) unusually hard and bright emission from 
one of the reddest Class I protostars in the c loud, Chandra sour ce 
#32 (MC09 source IC1396A-7; source e from lReach et al.N2004h '). 
Detailed X-ray spectral and variability analyses from IC1396A-7 
are out of the scope of this paper and will be presented elsewhere. 
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Figure 8. IR color-color diagrams of X-ray sources. For all panels, 
YSO diskfree (Class III) stars with available IR photometry in- 
formation are in red, YSO disky (Class I and II) stars are in blue, 
possible foreground stars unrelated to the region are in green. On 
panels (b)-(d) the redd ening vectors show A x ~ 2 mag using 
the extinction law from iFlahertv et"al | ||2007|) . (a) J - H versus 
H — K s . The dashe d line indicates the classical T Tauri star locus 
dMever et al.lll997l) . The solid line is a reddening vector marked 
at intervals of Ay = 10 mag. (b) K s — [3.6] versus [3.6] — [4.5]. (c) 
[3.6] - [4.5] versus [4.5] - [5.8]. (d) [3.6] - [4.5] versus [5.8] - [8.0]. 
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Figure 9. IR color-color diagram for non- Chandra sources with 
reliable photometry identified through the analysis of Spitzer- 
IRAC data in the current work (J3J. Classification criteria for 
selecting disky YSOs, [3.6] - [4.5] > 0.2 mag and [5.8] - [8.0] > 
0.2 mag, are shown by the dashed lines. Forty two non- Chandra 
YSOs additional to the disky stars from MC09 are identified 
(filled circles). The red dening vector shows Aj^ ~ 2 mag using 
the extinction law from iFlahertv et al.l ll2007f ). 

These are tabulated in column 8 of Table As in the case 
of th e X-ray selected YSO s in the Cepheus star forming re- 
gion jGetman et al.ll2009l . their Figure 6), the distribution of 
the spectral index for IC 1396A X-ray YSOs has a bimodal 
shape with peaks at ao ~ —1.2 and —2.7 corresponding to 
diskbearing and diskless star samples, respectively (graph 
is not shown). In agreement with the SED-based classifi- 
cation, all but one diskless X-ray stars have SED spectral 
indices of ao < —2.1, and all disky X-ray stars have indices 
of ao > —2.1. 

In Figure [8] we compare the SED-based classifications 
of X-ray stars with the expected loci of YSO stars in 
NIR and MIR color - color diagrams (e.g. lAllen et alj 12004 
iMegeath etafl 12004 lHartmann et al.ll2005l ). The positions 
of stars in the color-color diagrams are in good agreement 
with our previous classification of stars — simple color crite- 
ria that disky stars have [3.6] - [4.5] > 0.2, [4.5] - [5.8] > 0.2, 
or [5.8] - [8.0] > 0.2, select nearly all of the SED-based clas- 
sified disky stars. A couple of stars with the 3.6 — 4.5 color 
near zero and with an excess at longer wavelengths could 
be transition disk objects, candidates to have inner holes in 
their disks rather than being diskless. 



4 NON- CHANDRA INFRARED-EXCESS 
MEMBERS 

For a given mass, the X-ray detection efficiency of Class II 
stars is somewhat lower than that of Class III stars 
l|Getman et al.ll2009l ; IsTelzer et al.ll201ll . and Figure [5] here). 
To compensate for this effect, this section presents addi- 



tional non- Chandra IR-excess selected disky stars. The MIR 
catalog from MC09 provides an excellent ( £|7.4[) selection of 
disky members near/inside the globule. The catalog also in- 
cludes all di s ky st ars n ear /inside the globul e repo rted by 
iReach et al.1 l|2004 ) and ISicilia-Aguilar et all l|2006h . How- 
ever, MC09 covers only one fifth of the ACIS field. Thus, we 
use the results from our Spitzer-IRAC data analysis ( i]2.2l) 
to produce an additional sample of non- Chandra disky stars 
for the remaining ACIS field and beyond. 

Our IRAC data analysis finds 292 sources not detected 
by Chandra and not published by MC09, with photometric 
errors < 0.1 mag in all four IRAC bands. Figure shows 
the [3.6] — [4.5] versus [5.8] — [8.0] diagram for these sources. 
Using the IR color-color classification criteria presented in 
H3.3I we identify 42 (out of 292) IR-excess stellar candidates. 
Accor ding to the classification scheme of iGutermuth et al.l 
(2009), only one of these sources, source # 29, could be 
an extragalactic contaminant. Since the IC 1396 region is 
in the outer Galaxy, we do not expect the sample to be 
contaminated by asymptotic giant branch or classical Be 
stars. 

Table 2]reports their NIR and MIR magnitudes, appar- 
ent SED slope from IRAC photometry, a flag indicating if 
a source lies within the ACIS field, and stellar counterparts 
from SA05 and Bll. Out of 42 stars, 20 (without SA05 and 
Bll counterparts) are newly discovered members with 15 
of them found within the ACIS field. IR SEDs for these 15 
sources are given in Figure 1101 They deviate from photo- 
spheric SEDs at wavelengths longer than at least 3.6 /im. 
As expected, the majority of these new non- Chandra disky 
YSOs are fainter than the X-ray selected disky sample (com- 
pare cyan Ds and blue xs in Figure [6}. 



5 STELLAR AGE ESTIMATES 

Indications of an age gradient among PMS stars, decreas- 
ing from the central Tr 37 cluster towards the IC 1396A 
globule, have been reported by SA05 (see their Figure 11) 
and Bll (their Figure 17). We seek here to confirm this 
trend and find out if it is statistically significant by taking 
advantage of the increased stellar sample size near/around 
the globule and the availability of multi-epoch optical data 
for many sample members. This requires deriving and com- 
paring in a homogeneous way stellar ages of the combined 
optical SA05/SA10 and X-ray selected stellar sample. 

It is well known that estimates of individual stellar ages 
can be very uncertain due to several possible causes such as 
measurement uncertainty, stellar variability, binarity, differ- 
ent accretional histories, extinction unc ertainty, model un- 
certainty, and distance uncertain ty (e.g. iBaraffe et al. I l2009l ; 
ISoderblorn|[2010l ; |Preibisch|[20r2] ) . However, here we are not 
so much interested in precise ages for individual objects 
as in the general trends of age distributions and consis- 
tent treatment for comparison between different stellar pop- 
ulations. In spite of all the uncertainties, the methods of 
isochrones and color-magnitude diagrams often allow dis- 
criminatio n among YSO regions of different charact eristic 
ages (e.g. lBurningham et al.[|2005l ; iMavne et al.ll2007l ). The 
isochronal ages were shown to be in agreement with gravity 
related indices l|Lawson et al.|[2009l ). Similar isochronal ages 
found for binary components in the Taurus- Auriga region 
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Table 4. Non-Chandra disky stars additional to MC09 stars. Column 1: Source number for non- Chandra IR-excess sources found in our 
IRAC data analysis. These sources are located around the globule and are complementary to the IR-excess sources in/near the globule 
cataloged by MC09. Columns 2 and 3: IRAC right ascension and declination for epoch J2000.0 in degrees. Columns 4-7: 2MASS JHK S 
magnitudes, and 2MASS photometry quality and confusion-contamination flag. Columns 8-12: IRAC magnitudes and the aperture flag 
derived in this work. The flag gives photometric apertures and level of source contamination from nearby sources and nebular IR emission: 
2-pixel (1.73"), 3-pixel (2.6"), and 4-pixel (3.46") apertures with contaminating flux from a neighboring source of no more than 10%, 5%, 
and 5% in [3.6] band, respectively. Column 13: Apparent SED slope from IRAC photometry with ltr error. Column 14: A positional flag 
indicating if a source lies within the ACIS field of view ("1") or not ("0"). Columns 15 and 16: Stellar counterparts from Bll and SA05. 
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880500 
887000 
243333 
244708 
246292 
283917 
300708 
307917 
309125 
314375 
321458 
333417 
339625 
350333 
351917 
364167 
370583 
392083 
398792 
418667 
428667 
428792 
442167 
445292 
449833 
452208 
452292 
453917 
454958 
455958 
463292 
468958 
496833 
514542 
535125 
535625 
538667 
542417 
550000 
594958 
614250 
616833 



Dccl. 
(deg) 
(3) 



524333 

533194 1 

650250 

646667 1 

651556 

604472 

457306 

457528 

604889 

454722 

479806 

492861 

436778 

403167 1 

526639 

523222 

601222 

575333 

549556 

575889 

579500 

608750 

574444 1 

532278 

545083 

499667 1 

405056 

389111 

464278 

564333 

410000 

432250 

604472 

693028 

446583 

618778 

557250 

452333 

416861 

671500 

518917 1 

512889 



.7 

(mag) 
(4) 



14.53 ± . 
4.64 ± 
6.48 ± 

4.14 ± 

3.57 ± 

4.58 ± 
3.39 ± 
4.27 ± 
4.58 ± 
4.46 ± 
5.02 ± 
4.88 ± 
4.23 ± 

4.20 ± 
4.54 ± 

4.38 ± 
3.99 ± 
5.00 ± 
4.86 ± 
4.99 ± 
4.96 ± 
0.57 ± 
4.92 ± 
4.85 ± 

4.15 ± 
5.04 ± 

4.66 ± 0.1 

6.04 ± 0. 

6.39 ± 0. 
5.48 ± 0.1 

5.21 ± 0.1 
3.32 ± 0.1 
2.50 ± 0.1 
0.51 ± 0.1 
4.31 ± 0.1 
2.39 ± 0.1 

5.05 ± 0.1 
4.53 ± 0.1 
4.83 ± 0.1 
2.36 ± 0.1 
2.91 ± 0.1 



04 13 

03 12 

04 13 

03 13 

04 13 

05 14 
04 14 
07 13 



04 13 

04 14 

05 14 
04 14 



(mag) 
(5) 



63 ± 0.05 
81 ± 0.04 
40 ± 0.11 
34 ± 0.04 
65 ± 0.03 
65 ± 0.04 
21 ± 0.03 
15 ± 0.04 
59 ± 0.04 
48 ± 0.04 
20 ± 0.05 
11 ± 0.06 
36 ± 0.07 
29 ± 0.08 
79 ± 0.05 
34 ± 0.04 
03 ± 0.03 
98 ± 0.04 

02 ± 0.05 

09 ± 0.06 
18 ± 0.05 

10 ± 0.03 

03 ± 0.05 
95 ± 0.04 
14 ± 0.03 
31 ± 0.06 



13.68 ± 0.06 
15.23 ± 0.12 
15.04 ± 0.10 
14.68 ± 0.10 

14.42 ± 0.07 
12.40 ± 0.03 
11.48 ± 0.03 
9.72 ± 0.03 

13.43 ± 0.03 

11.57 ± 0.03 
13.83 ± 

13.58 ± 0.04 
13.78 ± 0.04 
11.36 ± 0.04 
11.98 ± 0.04 



13.22 ± 0.( 

13.51 ± 0.( 

14.62 ± 0.( 

12.97 ± 0.( 

12.26 ± 0.( 

13.27 ± 0.( 

11.52 ± 0.( 
12.40 ± 0.( 
13.01 ± 0.( 
12.95 ± 0.( 
13.93 ± 0.( 
13.79 ± 0.( 
12.89 ± 0.( 
12.82 ± 0.( 

13.28 ± 0.( 

12.98 ± 0.( 

12.63 ± 0.( 
13.49 ± 0.( 

13.62 ± 0.( 
13.77 ± 0.( 

13.76 ± 0.( 
9.84 ± 0.0 

13.64 ± 0.( 

13.63 ± 0.( 

12.77 ± 0.( 
14.05 ± 0.( 

13.25 ± 0.( 
14.30 ± .. 
14.28 ± 0.( 

14.24 ± 0.( 
13.97 ± 0.( 

11.89 ± 0.( 
10.85 ± 0.( 
8.77 ±0.0 

12.90 ± 0.( 
11.20 ± 0.( 

13.63 ± .. 
13.32 ± 0.( 
13.35 ± 0.( 
10.77 ± 0.( 

11.25 ± 0.( 



13 



UAAOcO 
AAA000|l3 
BB A 
AAA000I12 
AAA00011 
AAA00012 
AAA00010 
AAAcOO 
AAA00011 
AAA00012 
AAA00013 
AAA000|l3 
AAA. 
AAA. 
AAA000|12 
AAAcOO 
AAA00011 
AAA00012 
AAAcOO 
AAAOcO 
AAAcOO 
AAA000 8 
AAA00012 
AAA00012 
AAA00012 
AAA00013 

AAAcOO 

BBUcOO 

BAA000 

AAA00013 

AAA00013 

AAA00011 

AAA000 9 

AAA000 7 

AAA00011 
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AAA00012 

AAA00010 

AAA00010 



[3.6] 

(mag) 
(S) 



41 ± 
.00 ± 0. 
.79 ± 
.46 ± 0. 
.60 ± 
.60 ± 
.64 ± 0. 
.42 ± 
-91 ± 
.14 ± 
-31 ± 
.14 ± 
40 ± 
.94 ± 0. 
.59 ± 
.27 ± 
.86 ± 
.64 ± 0. 
.21 ± 
-13 ± 
-15 ± 
95 ± 0. 
.92 ± 
± 
.03 ± 
-19 ± 
.68 ± 
-37 ± 
.26 ± 0. 
.79 ± 
-34 ± 
.72 ± 
-13 ± 
83 ± 0. 
30 ± 0. 
.45 ± 
.60 ± 0. 
-10 ± 
.66 ± 
.36 ± 0. 
02 ± 
19 ± 



[4.5] 
(mag) 
O) 



12.15 ± 
12.64 ± 

13.23 ± 
12.21 ± 
11.12 ± 

12.21 Hz 
10.25 ± 
10.92 ± 

11.38 ± 

11.78 ± 
12.96 ± 

88 ± 
11.91 ± 
1.41 ± 

12.22 ± 
12.05 ± 
11.47 ± 

12.24 ± 
12.81 ± 
12.61 ± 

12.79 ± 
35 ± 0. 

12.56 ± 
12.55 ± 
11.69 ± 
12.78 ± 
13-08 ± 
11.98 ± 
13.90 ± 
13.24 ± 
12.87 ± 
13.18 ± 
10.73 ± 

9.51 ± 0. 

6.80 ± 0. 
11.01 ± 
10.22 ± 

12.39 ± 

12.40 ± 
11.98 ± 

9.65 ± 0. 

9.80 ± 0. 



[5.8] 
(mag) 
(10) 



12.06 ± 0.1 
12.16 ± 0.( 
12.66 ± 0-( 

11.90 ± 0.( 

10.57 ± 0-( 
11.84 ± 0.1 
9.96 ± 0.0 
10.52 ± 0-( 
10.95 ± 0.1 
11.42 ± 0-( 

12.55 ± 0-( 

12.58 ± 0-( 

11.57 ± 0-( 
10.84 ± 0.1 

88 ± 0-( 

11.58 ± 0-( 

11.09 ± 0.1 

11.91 ± 0.( 
12.45 ± 0-( 
11.95 ± 0.1 

12.32 ± 0-( 
7.59 ± 0.0 

12.20 ± 0.1 

12.21 ± 0.( 

11.36 ± 0-( 

12.37 ± 0-( 

12.59 ± 0-( 
11.64 ± 0.1 

13.33 ± 0-( 
12.57 ± 0.1 
12.37 ± 0.1 

12.56 ± 0-( 

10.36 ± 0.1 
9.24 ± 0.0 
6.41 ± 0.0 

10.78 ± 0-( 
9.85 ± 0.0 

12.30 ± 0-( 

12.10 ± 0.( 
11.74 ± 0.1 
9.33 ± 0.0 
9.62 ± 0.0 



[8.0] 
(mag) 
(11) 



11.73 ± 
11.11 ± 
11.57 ± 
11.24 ± 

9.48 ± 0. 
11.05 ± 

9.37 ± 0. 

9.78 ± 0. 

9.93 ± 0. 
10.86 ± 
11.55 ± 
11.94 ± 
10.43 ± 

9.69 ± 0. 
L0.87 ± 
11.15 ± 
10.43 ± 
11.45 ± 
11.68 ± 

10.47 ± 
11.62 ± 

15 ± 0. 
11.50 ± 
11.35 ± 
10.49 ± 
11.71 ± 
11.82 ± 
10.93 ± 
11.45 ± 

11.48 ± 

11.49 ± 
11.38 ± 

9.52 ± 0. 

8.44 ± 0. 

5.75 ± 0. 

10.18 ± 
3.95 ± 0. 
1.61 ± 

11.48 ± 

11.19 ± 
8.73 ± 0. 
8.81 ± 0. 



-2.04 ± 0.02 
-1.09 ± 0.02 
-0.27 ± 0.05 
-1.71 ± 0.01 
-0.84 ± 0.01 
-1.30 ± 0.02 
-1.42 ± 0.01 
-1.01 ± 0.02 
-0.79 ± 0.01 
-1.44 ± 0.02 
-0.93 ± 0.03 
-1.66 ± 0.03 
-0.51 ± 0.04 
-0.12 ± 0.05 
-1.22 ± 0.02 
-1.52 ± 0.08 
-1.33 ± 0.01 
-1.44 ± 0.02 
-1.11 ± 0.03 
0.34 ± 0.04 
-1.22 ± 0.03 
-0.27 ± 0.01 
-1.30 ± 0.02 
-1.36 ± 0.02 
-1.47 ± 0.02 
-1.16 ± 0.03 
-0.67 ± 0.02 
-1.28 ± 0.02 
0.68 ± 0.04 
-0.17 ± 0.05 
-0.90 ± 0.02 
-0.23 ± 0.03 
-1.21 ± 0.01 
-1.59 ± 0.01 
-1.03 ± 0.01 
-1.37 ± 0.01 
-1.15 ± 0.01 
-1.49 ± 0.04 
-1.71 ± 0.03 
-1.52 ± 0.02 
-1.44 ± 0.01 
-1.32 ± 0.03 
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are in agreemen t with the commonly assum ed coevality of 
binary systems l)Kraus fe Hillenbrand! 120091 ). 

First, we clean our source sample to reduce some of 
the effects that cause age uncertainty. Our initial sample is 
limited to the stars that have both FLWO and LAICA pho- 
tometric measurements. These are 89 X-ray stars and 53 
SA05-SA10 optical stars without X-ray counterparts. Fig- 
ure [11] (a),(b) compares their l/J-band magnitudes. The 
systematic difference between the Ic and Zj-band magni- 
tudes is due to the different photometric systems, FLWO- 
Cousins versus LAICA- Johnson. The linear transformation 
Ic = 0.2954 + 1.0138 x Ij minimizes the differences for 
all but a few sources with Ij > 17 mag. After this I- 
band correction, large spreads in AV and AI measure- 
ments are seen. The spreads can be mainly attributed to 
observational uncertainty and stellar variability. To suppress 
these effects, the sample is truncated to remove sources with 
|AV| > 0.4 mag and |A7| > 0.4 mag. The sample is further 
reduced by removing several sources with possibly high ex- 
tinction (J — H > 1.2 mag corresponding to Av > 3 mag) 
and/or with known spectral types G-Kl. For the latter, 
their measured ages could be ov erestimated by in correct con- 
sideration of birth line effects l)Hartmannl 12003). Our final 
source sample comprises 102 stars (black points in Figure 
HD(a),(b)). 



The apparent V-band magnitudes and V — Ic colors 
(and transformed V — /./ colors) are then de-reddened by an 
average extinction of < Av >= 1.56 mag (SA05,B11) usin g 
the reddening law of A^Av = 0.6 l|Bessell fc BrettJ Il988ll 
Comp arison with the PMS evolutionary tracks of lSiess et al.l 
|2000l ) yields the individual age estimates (tFLWo, thAiCA) 
shown in Table [3] 

For the final source sample, statistical photometric 
uncertainties for both V and /-band magnitudes are < 
0.04 mag. However, the systematic uncertainty still prevails 
- only 9% (22%) of the sources have AV < 0.04 mag 
(AI < 0.04 mag), but 68% (90%) have AV < 0.2 mag 
(AI < 0.2 mag). The FLWO-LAICA age differences quan- 
tify the effect of variability on the inferred age estimates: 
\tFLwo — tLAicA\/tFLwo is better than 10%, 20%, 30%, 
40%, 50%, and 60% for 27%, 49%, 60%, 70%, 75%, and 80% 
of the sources, respectively. 

Figure[TT](c) shows an optical color- magnitude diagram 
for the FLWO data (the LAICA data look similar). There 
are a dozen diskless stars with inferred ages > 10 Myr. Some 
of them could be members of a possible old distributed pop- 
ulation. Half have inferred ages exceeding 100 Myr. This 
value seems to be unrealistically high considering that the 
average age for young stars in the region is around 4 Myr 
(SA05). This may point to variability or to misclassified fore- 
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Figure 10. IR SEDs for 15 new non- Chandra disky YSOs within the ACIS field that are not in the source catalogs of SA05, MC09, or 
Bll. See Figure [7] for details. 



ground and background stars. For the analysis below, we 
limit the sample to stars with a more physically realistic 
age of < 10 Myr. 

Figure [12] shows the spatial distribution of t < 10 Myr 
stars stratified by the ages obtained from the FLWO data; 
the LAICA data produce similar strata. The region is 
further divided into two sub-regions: the East sub-region 
near/around the 06 ionizing star HD 206267 is associated 
with the central stellar cluster Tr 37, while the West sub- 



region near/around the IC 1396A globule is expected to 
be associated mainly with stars that have been recently 
formed in the globule. Our data provide strong evidence 
that the stars in the East sub-region are systematically 
older than the stars in the West sub-region (Figure [H] (d)). 
Kolmogorov-Smirnov (K-S) and "Fisher's exact" tests give 
small probabilities that the two age samples are drawn from 
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Figure 11. Derivation of ages. Comparison of V-band (a) and /-band (b) magnitudes for young SA05/SA10 and X-ray stars that have 
both FLWO and LAICA photometry. Black (red) circles mark 102 (40) stars that have V, /-band magnitudes consistent (inconsistent) 
between the FLWO and LAICA observations. Unity lines are shown in green, (c) VI C color-magnitude diagram for the 102 stars with 
consistent photometry: 61 X-ray selected stars (x) and 41 SA05/SA10 stars without X-ray counterparts (o). Colors indicate different 
disk classes: diskfree (red), diskbearing (blue), transition disk objects (green), and uncertain class (black). Magnitudes and colors have 
been corrected for the average extinctio n of < Ay >~ 1.5 6 mag (SA05). The black solid lines show PMS isoch rones for ages 1, 2, 3, 4, 
10, 100 Myr from evolutionary tracks of S iess et al. (2000) assuming distance of 870 pc llContreras et al.ll2002T) , The black dashed line 
shows a reddening vector of Ay = 1.56 mag. 



the same distributioijf]. The derived mean/median FLWO 
(LAICA) ages are 3.9 Myr/3.6 Myr and 2.8 Myr/2.2 Myr 
(4.4 Myr/4 Myr and 2.7 Myr/1.9 Myr) for the East and 
West sub-regions, respectively. These are consistent with the 
previously reported average age of 4 Myr for the Tr 37 clus- 
ter (SA05). Typical age estimates for the stars in front of 
the globule are t < 2 — 3 Myr. Meanwhile, younger ages 
(t < 1 Myr) should be assumed for numerous disky objects 
that are embedded in the globule and lack optical counter- 
parts (JO}. 



6 SPATIAL STRUCTURE 

In this section we combine data from different stellar cata- 
logs in order to examine spatial distributions of all known 
young stars in/around the globule and the neighboring Tr 37 
cluster. For the central part of the Tr 37 cluster, we uti- 
lize optical data from SA05, while for region in and around 
the globule we consider the MIR catalog of MC09, the X- 
ray catalog, and our additional MIR catalog from Table [4] 
The latter t h ree ca talo gs in turn include disky YS Os from 
iReach et al.1 i|2004! ) and ISicilia-Aguilar et all (|2006l ) as well 
as all Ha stars recently identified within the ACIS field by 
Bll. See Figure [T3l for the spatial distribution of these stars. 



6 In the case of the FLWO ages, P KS = 1.6% and a 2 X 2 
contingency table with younger/older stars in the East: West 
region = 14/16 : 42/15 is consistent with PFisher = 1-8%. 
In the case of the LAICA ages, P K s = 0.01% and a 2 X 2 
contingency table with younger/older stars in the East: West 
region = 11/18 : 45/14 is consistent with PFiaher — 
0.08%. The derived probabilities are significantly lower than 
the conventional criterion for statistical significance of < 5% 
(http:/ /en. wikipedia.org/wiki/Statistical_significance ). Thus the 
null hypothesis that the two samples are drawn from the same 
age distribution can be confidently rejected. 



6.1 Stellar surface density gradient from the 
SA05 sample 

Here we define the central portion of the Tr 37 cluster as a 
circular area of radius i? = 15' (large red circle in Figure [T5J 
centered on the 06 star HD 206267 and bordering the head 
of the IC 1396A globule. We refer to this area hereafter as 
R15. Nearly the entire area is within the central 4-Shooter 
field (Figures 1, 9 in SA05), which has a uniform and deep 
optical exposure. The circular area comprises 150 optical 
young stars (Tables 1, 3, 5 in SA05). The cluster is not cen- 
trally concentrated. For instance, compared to the younger 
Orion Nebula Cluster (ONC) where the surface stellar den- 
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Figure 12. Age analysis for the combined optical SA05 and X-ray selected young stellar sample. Age-stratified stellar sub-samples are 
shown on different panels: t < 2 Myr (a), 2 < t < 3.5 Myr (b), and 3.5 < t < 10 Myr (c). The polygon indicates the ACIS field of view. 
The dashed line at roughly half projected distance between the 06 star HD 206267 and the IC 1396A globule is a demarcation line that 
divides the region into two sub-regions, East and West, (d) Cumulative distributions of age are compared between the 30 East and 57 
West stars with available ages. 



sity dr ops by a factor of > 25 from the center to the 2.5 pc 
radius (|Hillenbrand fc Hartma nn 1998), the density in Tr 37 
decreases by only a factor of two over a similar parsec scale. 

Our goal here is to search for local stellar surface over- 
densities on the periphery of R15 area that could be re- 
lated to triggered star formation in now-dissipated molecu- 
lar cloudlets. The overlap region between R15 and the ACIS 
field (labeled as "C5" in Figure [T3] (a)) is especially useful 
in regards to possible triggered star formation in the origi- 
nally larger IC 1396A cloud. The periphery of R15 is divided 
by five control equal-area segments labeled as "C1"-"C5". 
From the SA05 master sample, we count 16, 20, 18, 15, and 
35 stars in segments CI, C2, C3, C4, and C5, respectively. 
We thus find that the stellar surface density in C5 area is 
two times higher (with statistical significance P > 99.5%) 
than the average density in the remaining periphery area of 



the central cluster. The ~ 17 extra stars could be associ- 
ated with triggered star formation in the originally larger 
IC 1396A cloud. 



6.2 Increase in the census of YSOs in and around 
IC 1396A 

As X-ray emission from PMS stars is based on enhanced 
solar-type magnetic reconnection events rather than disk or 
accretion processes, X-ray selection delivers rich and clean 
samp les of diskless stars missed by IR selection l|Feigelsonl 
2010). Previous studies have identified 129 YSO members 
within our ACIS field of view. This study has doubled that 
census by identifying 130 additional stellar members in our 
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Figure 13. Spatial distribution of known young members of the region over-plotted on DSS2-Red (a) and Spitzer-IRAC 8 fim (b-d) 
images. On all panels the big red circle outlines the area of the central cluster Tr 37 of radius R = 15' centered on the 06 star HD 206267. 
The blue polygon outlines the ACIS field. The green polygon demarcates the globule, tracing 8.0 /an emission from hot dust and PAHs. 
(a) Positions of young members from the SA05 (red o) and the X-ray (blue x) catalogs. "CI" through "C5" label five equal-area segments 
at the edge of the central cluster. Positions of the two optically bright rims in front of the cloud are marked by arrows and labeled in 
green, (b) Positions of all known young members of the region within the ACIS field: optical SA05 (red o), X-ray selected (blue x), MIR 
stars from MC09 (green □), H« stars from Bll (magenta +), and additional MIR selected stars from this study (cyan o). (c) Separation 
on X-ray selected diskless stars (red x) and disky stars (blue +) combined from MC09 and this work, (d) Contour maps of source surface 
densities from panel (c) smoothed with a Gaussian a = 1.4'. Source density contours are in linear units from 0.75 to 1.5 sources per 
arcmin 2 for diskless (red), and from 0.75 to 3 sources per arcmin 2 for disky (blue) stars; both with a step of 0.25 sources per arcmin 2 . 
Stellar enhancements above 1.25 source per arcmin 2 level are indicated by arrows and labeled in blue. Known YSOs with spectral type 
earlier than B6 are marked in cyan. North is up, east is to the left. 



X-ray ( %T2t and MIR (© catalogs (Figure[H(b)). Twenty- 
nine (14 Chandra and 15 non- Chandra) exhibit an IR excess. 



6.3 Stellar surface density gradient from the 
combined stellar sample 

A map of the combined sample of 259 members is shown in 
Figure [13] (c) . They appear clustered rather than uniformly 
distributed, and the distributions of the disky (blue +) and 
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diskless (red x) stars appear to be offset horizontally from 
each other. This is more plainly seen in the source surface 
density maps shown in panel (d), which were constructed 
using a fixed Gaussian kernel. 

Directly in front of the molecular globule, most of the 
diskless and a portion of the disky stars lie in a > 3 pc long 
and ~ 0.7 pc wide layer. This result is robust against ob- 
servation selection effects as the excess of the stars within 
the layer is seen for both the Chandra X-ray selected disk- 
free stars and the disky stars that were mainly identified 
through Spitzer IR selection. The position of the layer is 
also consistent with the positions of the very young optical 
stars (SA05, Bll) and with the excess of the optical stars 
on the western border of Tr 37 that was discussed in £16. II 

The layer can be further subdivided into a few stellar 
groupings: about 30 disky and diskless stars along the bright 
rim and to the east of the globule towards the direction of the 
major ionizing system HD 206267 (labeled as "Gl"); about 
20 disky and diskless stars to the north-east of the globule 
around the B3+B5 high-mass system HD 206081 ("G2"); 
about 20 stars to the south-east of the globule (the disky 
"G6" and the diskless "G7"). Some groupings are disky and 
others diskless, suggesting different ages from distinct RDI 
events in different cloud cores over millions of years ( H8.3.2p . 

Inside the globule, the stellar population is dominated 
by disky objects, whose spatial distribution is also highly 
non-uniform and clumpy, with the bulk of the stars lying 
along the western edge of the globule that faces the ionizing 
system HD 206267. Two major stellar groupings are seen 
inside the cloud: a larger group composed of ~ 50 Class 
III/II/I sources at the head of the globule behind the bright 
optical Rim A (group "G3"); and a smaller group composed 
of ~ 10 sources, mainly Class I protostars, located behind 
optical Rim Aa in the portion of the cloud named SFO 36 
(group "G4"). 

6.4 Disk-fraction gradient 

The source surface densities of the diskless and disky stars 
lying within the ACIS field inside and in front of the cloud 
suggest a spatial gradient of apparent disk fraction (Figure 
Q3](d)). The highest disk fractions are inside the cloud, be- 
hind Rims A (9/9 = 100%) and Aa (35/47 = 75±6%f|. An 
intermediate disk fraction of 41/85 = 48 ± 5% is found in 
the layer with the enhanced stellar surface density in front of 
the globule (groups "Gl", "G2", "G6", and "G7" =G1267). 
Eastward of the layer, towards Tr 37, the disk fraction is 
lowest (12/34 = 35±8%) and is consiste nt with a recent est i- 
mate of the intrinsic disk fraction (39% lM~ercer et al]|2009|) . 
This disk-fraction gradient within our sample is a tenta- 
tive finding, for two reasons. First, the gradient is statis- 
tically significant across the G1267-globule region, but not 
across the Tr37-G1267 regiorjf]. Second, the comparison is 
performed with apparent disk fractions, which are subject 

7 Errors on disk fractions have bee n estimated using bino mial 
distribution statistics as described by [gurgasscr et al.l II2003T ) . 

8 For the G1267-globule area, a 2 X 2 contingency table with 
disky /diskfree ratio = 41/44 : 35/12 is inconsistent with the null 
hypothesis of no disk-fraction gradient (PFisher — 0.5%, using 
Fisher's exact test). However, for the Tr37-G1267 area, a 2 X 2 
contingency table with disky /diskfree ratio = 12/22 : 41/44 is 



to a spatially varying X-ray sensitivity that is suppressed 
inside the globule. 



7 XLFS, IMFS FOR THE STELLAR 

POPULATIONS IN THE TR 37/IC 1396A 
REGION 

We construct X-ray luminosity functions (XLFs) and ini- 
tial mass functions (IMFs) independently for the different 
projected regions on the sky. The separation into the older 
Tr 37 and younger (likely triggered) IC 1396A stellar popu- 
lations will be made in iJS] For stars outside the globule in 
the ACIS field, two independent analysis are performed: an 
XLF analysis for purely X-ray selected stars ( JZ7TJ , and an 
IMF analysis for the combined X-ray and MIR YSO sample 
(SJ73J. The IMF analysis of the optical SA05 YSO sample 
is employed to measure the total stellar population within 
the R = 15' radius circular area around HD 206267 (E3)- 
Finally, the IMF of a stellar sample mainly composed of 
disky YSOs from MC09 combined with some additional X- 
ray stars is utilized to estimate the total number of stars 
inside the globule within the ACIS field ( 872)1 ■ 



7.1 XLF of the stellar population around the 
globule 

Comparison of X-ray luminosity functions (XLFs) between 
newly studied young stellar clusters like IC 1396A and 
a well-characterized cluster like the ONC can give in- 
sight into any missing p opulation of stars. The ass ump- 
tion of a universal XLF (|Feigelson fc Getman I I2005T ) has 
been used in the past to estimate tot al populations in sev - 



eral y oung clusters, suc h as Cep B (iGetman et al.l 



NGC 2244 (Wang et al 



M17 llBroos et all I2007TI NGC 6357 JWang et al 



l2010l ). W40 (Kuhnetal 



2006), 



2007), 



20081) . NGC 2237 (I Wang et al.1 



20 id ), Tr 15 l|Wang et al.1 1201 ll ). 



Tr 16 (IWolketal.1 120111 ). Sh 2-254/ 255/256/257/258 
jMucciarelli et al.ll201 ll ). and NGC 1893 (jCaramazza et all 
120111 ). Notice that XLFs are relatively insensitive to age ef- 
fects; between 0.1 and 10 Myr, only a slight change (roughly 
0.3 in logLx) is found in the ONC stellar po pulation 
|Preibisch fc Feigelsonll2005l : iPrisinzano et al.ll2008l ). 

For the X-ray selected young stars related to IC 
1396A, X-ray luminosities (Ljt ) and absorbing column den- 
sities (N h ) are generated usin g the non-parametric method 
XPHOT jGetman et alj|2010l ). The concept of the method 
is similar to the long-standing use of color-magnitude dia- 
grams in optical and infrared astronomy, with X-ray median 
energy replacing color index and X-ray source counts re- 
placing magnitude. Empirical X-ray spectral templates de- 
rived from bright sou rces from the Chan dra Orion Ultra- 
deep Project (COUP; IGetman et ai1l2005h are further used 
to translate apparent photometric properties of weak PMS 
stars into their intrinsic properties. The advantage of the 
XPHOT method over a traditional parametric spectral mod- 
eling (e.g., XSPEC) is that it is more accurate for very faint 
sources and provides both statistical and systematic (due to 



consistent with the null hypothesis of no disk-fraction gradient 

(PFisher = 23%. 
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Figure 14. Comparison between XLFs of the ONC and the pop- 
ulation of the X-ray selected YSOs located outside the IC 1396A 
globule. The dashed-dotted histogram shows the COUP unob- 
scured cool sample of 839 stars; the solid histogram is for the 
sample of 91 X-ray YSOs with available X-ray luminosity infor- 
mation located around the IC 1396A globule. Error bars indicate 
68% confidence intervals (1 <r) from Monte Carlo simulated distri- 
butions as described in the text. The thin dashed lines are added 
to aid the eye; they are based on the shape of the ONC XLF and 
scaled downward to match the IC 1396A XLF. 
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Figure 15. The mass distribution of all known X-ray and IR- 
excess selected YSOs outside the globule within the ACIS field 
of view, assuming an age of ~ 2 Myr. Red and blue histograms 
are for diskless and disky YSOs, respectively. The lcr Poisson 
errors (Gchrcls 1986) are also shown. The black line repr esents 
the scaled version of the Galactic field IMF l|Kroupall2002 ). 



tion do not affect the results of the analysis because their 
X-ray luminosities are expected to be below the complete- 
ness limit. 



uncertainty in X-ray model) errors on derived intrinsic fluxes 
and absorptions. The XPHOT results are given in Table [3] 

Figure [14] shows the XLF for all X-ray selected YSOs 
with available Lx information that lie outside the IC 1396A 
globule. The error bars on the XLF were generated from 
1000 Monte-Carlo simulated XLF distributions with indi- 
vidual X-ray luminosities randomly drawn from Gaussian 
distributions with mean equal to the measured source's 
\og(Lx) and variances based on the corresponding statistical 
and systematic errors summed in quadrature. Over-plotted 
is a template XLF fo r the COUP unobsc ured population 
(839 cool ONC stars in lFeigelson et al.ll2005h. which is nearl y 
complete down to M = 0.1 M (§7.2 in lGetman et alf 2 006V 

The IC 1396A XLF falls quickly below log(Lx) < 
30.0 erg s _1 due to incompleteness, with a conservative 
value for the completeness limit at arou nd log(Lx) = 
30.5 e r g s -1 . From the M-Lx correlation of IPreibisch et all 
(2005); iTelleschi et al.1 ((2007) this corresponds to a mass 
value of M ~ 1 Mgfl Matching the two cluster XLFs in 
the range 30.5 < log(Lx) < 31.5 erg s _1 requires scaling 
down the ONC unobscured population by a factor of ~ 0.28, 
yielding a total of ~ 235 stars down to 0.1 Mq outside the 
IC 1396A globule within the ACIS FOV. Notice that the 
observed weak X-ray sources without available Lx informa- 



9 The alternative IMF analysis f £|7.2| I will further show that the 
inclusion of additional non- Chandra Class II members selected by 
IR-excess will push the mass completeness limits for all known 
YSOs outside the globule down to 0.3 - 0.4 Mq. 



7.2 IMF of the stellar population around the 
globule 

In the area outside the globule covered by the ACIS field and 
all four bands of the IRAC mosaic, we examine the stellar 
IMF separately for diskless and disky stars. To compensate 
for the effect of higher X-ray detection efficiency of Class III 
stars, we use the disky stellar sample that is a combination 
of the Chandra and non- Chandra IR-excess member young 
stars (Tables |3] and g] in this work, and Table 2 from MC09) 
while retaining the diskless stellar sample as purely com- 
posed of X-ray stars (Table [3]). 

Approximate stellar mass estimates are obtained based 
on star locations in the J versus J — H color-magnitude 
diagram and theore tical stellar model tracks derived by 
Baraffe et ail dl99ct) (for 0.02 M© sC M < 1.4 Mq) and 
Siess et al.l l|200fj| ) (for 1.4M <M< 7.0 Mq). These color- 
magnitude diagrams are shown in Figures [5] and [6] We 
are aware that these mass estimates are subject to signif- 
icant uncertainties and may be incompatible with individ- 
ual masse s obtained by othe r methods such as optical spec- 
troscopy l|Kuhn et al.ll2O10l . their Appendix section). How- 
ever, here we are not so much interested in exact masses 
for individual objects as in the general trends of mass dis- 
tributions and consistent treatment for comparison between 
different stellar populations. 

Typical derived age estimates for the stars in front of 
the globule are t < 2-3 Myr (©. Figure [151 shows the IMF 
of YSOs outside the globule for an assumed age of 2 Myr. 
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Figure 16. The NIR color-magnitude diagram for all (inside and 
outside the ACIS field) optical YSOs from SA05. Low-mass stars 
are shown as black circles; intermediate and high-mass stars are 
shown as green triangles. The unabsorbed zero-age main sequence 
and 4 Myr PMS isochrone at 870 pc distance are shown as dashed 
blue and solid magenta lines, respectively. Reddening vectors of 
Ay = 5 mag for masses 0.5 and 1 Mq are indicated by dashed 
magenta lines. Stars of spectral type earlier than B6 are outlined 
by red □ and labeled. 



The mass distribution within errors follows the shape of the 
Galactic field IMF with rough mass completeness limits at 
0.3—0.4 Mq. The disk fraction obtained by counting diskless 
and diskbearing stars with masses above the completeness 
limit is 41 ±5%. The inferred total number of stars is around 
250 stars down to 0.1 Mq. Similar results are obtained from 
the IMF analysis with the assumption of an age of 3 Myr 
(figure is not shown). The total number of stars (~ 250) 
obtained here from the IMF is consistent with the total (~ 
235) derived from the age-insensitive XLF analysis ( 97,1( 1, 



7.3 IMF of the central part of Tr 37 

We estimate the total number of young stars in the central 
part of Tr 37 within the R = 15' circular area centered on 
the 06 star HD 206267 using the optical catalog of SA05. 
For the central cluster we assume the age of ~ 4 Myr (SA05, 
and ^S] here). As in 97.21 stellar mass estimates are derived 
from the J versus J — H color-magnitude diagram. 

Figure [16] shows that about 20 intermediate-mass YSO 
candidates lie far to the left from the 4 Myr isochrone with 
some of them along the ZAMS locus. It is unclear if these 
are members of an old distributed population or older stars 
unrelated to the region. These stars are excluded from the 
IMF analysis. 

Figure [T7I shows that the mass distribution for the opti- 
cal stars in the central cluster within errors follows the shape 
of the Galactic field IMF. The approximate mass complete- 
ness limit for the optical sample is 0.8 Mq. The inferred 
total number of stars in the circular R — 15' cluster area 
is found to be 480 stars down to 0.1 Mq. This IMF popu- 
lation and the XLF population from 97.11 are consistent in 
the overlap region between the circular R — 15' cluster area 
and the ACIS field. 
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Figure 17. The mass distribution for optical SA05 stars located 
within the circular area of radius R = 15' around the 06 star 
HD 206267 (solid histogram). The das hed line repre sents the 
scaled version of the Galacti c field IMF dKroupalbOol) . The la 
Poisson errors <Gehrelj|l986l) are also shown. 



7.4 Estimate of the stellar population inside the 
globule 

The J versus J — H diagram of X-ray selected YSOs shown 
in Figure [5] suggests that the X-ray detection efficiency of 
disky stars is lower than that of diskless stars. The J versus 
J — H diagram for all known disky YSOs with available NIR 
information (Figure [6} further shows that Chandra did not 
detect many low-mass and/or highly obscured disky stars 
identified by means of Spitzer data. For the analysis of the 
stellar population inside the globule, a region where Chan- 
dra's sensitivity deteriorates, we rely heavily on the Spitzer 
catalog of disky stars from MC09. 

There are 18 X-ray selected diskless stars and 45 IR- 
excess selected and/or X-ray selected diskbearing stars pro- 
jected against the globule (Table [3] here and Table 2 from 
MC09). With the assumption that all of these observed stars 
are globule members, the apparent disk fraction averaged 
across the entire globule is 45/(18 + 45) or 70%. Twenty of 
the diskbearing stars have been classified as Class I objects 
(MC09). The radiative transfer model fits to the optical-IR- 
radio band data of a mid-IR selected sample of YSOs in the 
globule yield ages of < 0.2 Myr an d < 1 Myr for the Class I 
and Class II objects, respectively ( Reach et al.| [2009). Con- 
sidering the high apparent disk fraction, the high fraction of 
Class I objects, and the young ages derived for th e selected 
sample of Class I /ll YSOs by iReach etHl < |2009h . we rea- 
sonably assume an age < 1 Myr for the entire population in 
the globule. 

The IMF analysis here is restricted to the diskbearing 
population in the globule. These are 38 out of 45 observed 
disky objects with available stellar mass estimates obtained 
from the J versus J — H color-magnitude diagram (Figure 
[6]). Their mass distribution is shown in Figure [TBI Despite 
the low counting statistics, the shape of the distribution is 
reminiscent of that of the Galactic field IMF. The IMF anal- 
ysis gives a crude number on total population of diskbearing 
stars in the globule of ~ 40 stars down to 0.1 Mq. This is 
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Figure 18. The mass distribution of disky stars found within the 
globule. Mass distribution of thirty eight young disky stars (36 
from MC09 and additional 2 new from the X-ray catalog) with 
available mass estimates, which lie projected against the globule, 
is shown as a solid histogram. The das hed line repre sents the 
scaled version of the Galacti c field IMF l lKroupal I200I) . The la 
Poisson errors jGehrelslll98ll are also shown. 



likely a low limit because 7 diskbearing stars without JHK 
information have not been included in the analysis. 

With respect to the diskless population in the globule, 
18 X-ray selected dikless stars projected against the globule 
have not been included in the IMF analysis. Out of these 18 
X-ray stars, at least 10 have relatively high source extinc- 
tions, comparable to those of the diskbearing population in 
the globule (Figure [SJ. Assuming an age of 1 Myr, all 10 
would have masses above 0.1 M@. Thus, the total number 
of disky and diskless stars in the part of the globule covered 
by the ACIS field of view is > 60. 



8 IMPLICATIONS FOR TRIGGERED STAR 
FORMATION 

8.1 Observational evidence for RDI triggered star 
formation 

Spatial-age gradients found in stellar populations in/around 
BRCs are often viewed as strong observational evidence 
for triggered star formation. The youngest stars are typ- 
ically embedded at the edge of the cloud facing the ion- 
izing star, and older stars are aligned within and in front 
of the bright rim towards the ionizing star. For example, 
spatio-temporal gradients and triggered star formation have 
been reported in BRCs and pilla rs on the edges of large 
Hn r egions in the Eagle Nebula dSugitani et al. I 2002T ). IC 
1848 llMatsuvanagi et alj|2006t lOgura et all 120071). IC 1396 
llGetman et ail 120071 ; Ikeda et ail |200S|; IChoudhurv et al] 
|2010| ), Sharpless 15 5 jGetman et ail I2009T ), Sharpless ~171 
and Sh arpless 296 ( Chauhan et al. 20091 ). and the Carina 
Nebula (|Smith et alll201Ch . 

In Sj5]and[6]we establish the spatio-temporal gradient of 
young stars from the center of the IC 1396A cloud toward 
the HD 206267 06.5f star. The identified stellar members 
are not uniformly distributed around the globule but rather 
are concentrated along the ionized rim, both outside (stel- 



lar groups "G1" + "G2" + "G6" + "G7" in Figure [i3)l and in- 
side ( "G3" + "G4" ) the cloud. The clustering of young stars 
in front of the cloud is also consistent with the previously 
reported arc-shaped layer of ~ 20 low-mass high accretion 
stars from Bll and/or a dozen of ~ 1 Myr optical stars from 
SA05. The stars located inside the cloud have typical ages 
of < 1 Myr f £|7.4[l . the stars found right in front of the cloud 
have typical ages of < 2 — 3 Myr (SJSJ, and the stars located 
further to the east, in the Tr 37 cluster, are older with ages 
of > 3.5 Myr (*j5j. This age gradient is also consistent with 
our tentative finding that the apparent disk fraction is de- 
creasing from 75 — 100% inside the cloud cores, to 48% in 
front of the cloud, and to 35% further to the east towards 
HD 206267 (©. 

The results from this work, combined with the results 
from the previous studies, thus provide strong evidence that 
the Tr 37/IC 1396A region possesses many of the observa- 
tional features of an RDI triggered star formation: 

(i) The presence of an exciting star (HD 206267) and 
a bright-rimmed cloud (IC 1396A) protruding inside the 
Hn region and surrounded by an ionized rim (Rim A and 
Rim Aa) facing the exciting star (Figure [T}. 

(ii ) The presence o f dense molecular cor es close to the 
rim (iPatel et ail [19951 , their Figure 8) and dWeikard et ail 
Il99rj . their Figure 14a). This is consistent with the scenario 
where the expanding ionization front from HD 206267 pene- 
trates the inter-core cloud region and drives shocks into the 
cores, compressing them sufficiently to trigger gravitational 
collapse. 

(iii) A rough equilibrium between the Hn pressure and 
the cloud pressure found at the interfaces between the Rim A 
and the head of the IC 1396A globule (|Reach et ail I2009T) 
and the Rim Aa and the SFO 36 cloud l|Morgan et ail 
120041) . This indicates that the cloud could have been re- 
cently shocked by a photoionization shock. I n addition, for 
the SFO 36 cloud. iMorgan et ail (|2009l . l2010l ) measure rela- 
tively high values of CO excitation temperature and turbu- 
lent velocity dispersion, which are consistent with those of 
other shocked cloud candidates. 

(iv) A clumpy spatial distribution of young stars both 
inside and in front of the cloud (Figure 1 13 p suggesting re- 
cent formation in distinct molecular cores. The proximity of 
the stars to the ionized rim immediately suggests a physi- 
cal relationship between star formation and the ionization 
front propagating through the globule, and can be explained 
within the framework of the RDI model. 

(v) The spatio-temporal gradient of young stars, both age 
and disk-gradient, oriented towards the exciting star (Fig- 
ures [HEM here, Figure 11 in SA05, and Figure 18a in Bll). 

8.2 Contribution of triggered star formation to 

the total population of the IC 1396 H 11 region 

8.2.1 Motivation 

The overall effect of triggered star formation on the global 
star formation ra te of the Galaxy is difficult to ascertain 
l|Elme grccn 201l]). A handful of studies attempted to esti- 
mate the effect. From an observational side, counting IRAS 
point sources embed ded in BRCs of numerous Hn regions, 
ISugitani et ail (|l99ll ) estimate the contribution of triggered 
star formation from expanding H 11 regions to be only ~ 
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5% by mass of the stars in the Galaxy. iDeharveng et all 
\20im find that 26% of the mid-IR shells from the Spitzer- 
GLIMPSE survey harbor UCHlls and/or methanol masers, 
likely indicating triggering of massive star formation. From 
comparison of positions of very young and massive stars 
identified in the RMS survey in respect to locations of nu- 
merous mid-IR s hells f rom the Spitzer-GLIMPSE survey, 
I Thompson et all (|201 lh estimate the fraction of Galactic 
massive stars triggered by expanding Hll regions to be 
14 - 30%. 

From a theoretical perspective, by turning on and off 
an irr adiation fee d back o n a turbulent 10 4 Mq molecular 
cloud, iDale et all (|2007bf) find that triggered star forma- 
tion might increase the total amount of star formation by 
as much as 30%. However, their SPH simulations show no 
observable characteristics for distinguishing triggered from 
spontaneously formed stars. Although it is often difficult to 
make such a distinction in observational studies of individ- 
ual H n regions, in £|8. II and 18.2.21 we argue that an RDI 
triggered stellar population can be identified in the case of 
the Tr 37/IC 1396 A region. Hence, the contribution of trig- 
gered star formation in the central part of the IC 1396 H 11 
region, and with some additional assumptions in the entire 
Hn region, can be estimated. 

8.2.2 Contribution in the central part of the Hll region 

In [JT] we establish the total populations of stars down to 
0.1 Mq for different areas of the region: ~ 480 stars within 
the R = 15' radius circular area around HD 206267; ~ 235 — 
250 stars in the area outside the globule covered by the 
ACIS field; and > 60 YSOs inside the globule within the 
ACIS field. These numbers are employed here to estimate 
the input of triggered star formation to the total population 
of the central part of the region. 

Based on the arguments presented in £18.11 it is reason- 
able to propose that most of the YSOs inside the IC 1396A 
cloud and the bulk of the stars in the stellar groups "Gl" + 
"G2" + "G6" + "G7" (=G1267) in front of the cloud (Fig- 
ure I13[) were formed through the RDI process. The notion 
is also consistent with our finding of a 2-fold increase in 
the surface density of the optical SA05 stars at the west- 
ern border of the central Tr 37 cluster ( - Assuming that 
the 2-fold increase is fully due to a triggered population, we 
obtain ~ 20% as a conservative estimate on the fraction of 
the old stellar members of Tr 37 within the G1267 group, 
which implies that ~ 50 out of the 60 X-ray YSOs detected 
in G1267 have a triggered origin. 

The total number of triggered stars in G1267 group in 
front of the BRC is estimated to be ~ 80 = 50 x (235/146) 
stars down to 0.1 Mq, where ~ 50 is the number of Chandra 
observed triggered candidates in G1267, 146 is the number 
of Chandra observed stars in the entire ACIS field outside 
the globule, and 235 is the inferred total number of stars 
in the entire ACIS field outside the globule. This number 
combined with the conservative (low limit) estimate on the 
total stellar population inside the globule (> 60 stars; £I7.4[1 
gives the total triggered population in IC 1396A of > 140 
stars down to 0.1 Mq. This is likely a lower limit as we do 
not include possible contribution of the sparse group of stars 
with a high apparent disk fraction seen right in front of the 
Rim Aa (Figure ITSb). 



In contrast, the total population of the central Tr 37 
cluster is 480-35/2/150 x 480 = 424 stars down to 0.1 M©, 
where 480 is an estimate on the total stellar population 
within the R = 15' radius circular area around HD 206267, 
35 and 150 are numbers of the observed optical SA05 stars 
located within the overlap area "C5" and the entire R = 15' 
circular area, respectively, and 2 accounts for the increase in 
the surface density of the stars within the overlap area due to 
the presence of the triggered population. The contribution 
of the triggered stellar population to the total population of 
the central part of the IC 1396 H 11 region to date is thus 
>25% (> 140/(424 + 140)). 

8.2.3 Contribution for the entire Hll region 

A rough estimate of the contribution of triggered star for- 
mation for the entire IC 1396 Hll region can be further 
obtained by extrapolating our results for the triggered pop- 
ulation in/around IC 1396 A cloud combined with the infor- 
mation on the low-mass high-accretion stars from the survey 
of the entire IC 1396 Hll region by Bll. This requires the 
following additional assumptions. First, we assume that the 
total number of stars in the entire IC 1396 H 11 region is sim- 
ilar to that of the Orion Nebula Cluster. This is supported 
by the fact that they each have equal numbers of stars B9 or 
earlier, and thus should have roughly the same total stellar 
population^]- Here we assume that the ON C population is 
2800 stars l|Hillenbrand fc Hartmannlll998l ). Second, we as- 
sume that the birth of all low-mass high-accretion Bll stars, 
which are found in clusters in front of the BRCs, was trig- 
gered. This is true for the ACIS field, where out of 22 Bll 
stars 21 are likely associated with the triggered population: 
19 belong to the G1267 stellar group in/around Rim A, one 
is within the cloud's head, and the other one is right in front 
of the Rim Aa. Based on Figure 9 from Bll we count > 15 
YSOs clustered in front of Rim E, about 10 stars in front 
of Rim B, several high-accretion stars in front of Rims J 
and H, and > 10 stars in front of some other, unlabeled, 
BRCs, totaling > 40 observed Bll triggered stars outside 
the ACIS field. Third, we assume that the ratio of the num- 
ber of the observed Bll triggered stars to the total number 
of the associated triggered stars, Rt, is constant across the 
Hn region. Within the ACIS field, R t = 21/(> 140) ~ 15%. 
The contribution of triggered star formation for the entire 
Hll region is then > 14% ((> 40 + 21)/7?t/2800). 

Combining this with the estimate for the central cluster 
in 38^21 > 14 - 25% of the stars in the IC 1396 region 
have a triggered origin. This is consistent with the various 
observational and theoretical estimates reviewed in £|8 . 2 . 1 1 

However, it has to be stressed that we are observing 
IC 1396 and other giant Hll regions at only a single mo- 
ment in their histories. IC 1396 has expanded to a radius 
12 parsecs over few million years, and the Hll shock likely 

10 From Table 1 of lStelzer et al.l l l2005h . the ONC has 9 stars be- 
tween B3-07 and 5 stars between B9-B4. From Table 5 of SA05 
and adding HD 206267 and HD 206183, IC 1396 has 3 probable 
members between B3-06 and 13 between B9-B4. So they each 
have ~ 15 stars B9 or earlier. Furthemore, the most massive sys- 
tems in both regions have similar spectral types — 06 star HD 
206267 and 09.5 star HD 206183 in IC 1396 vs. 06-7 star 6 1 Ori C 
and 09.5 star 6» 2 Ori A in ONC. 
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triggered star formation in many now-dissipated molecular 
cloud cores. The H II region is likely to persist for several mil- 
lion years into the future with supernova remnant shocks in 
addition to ultraviolet radiation and OB wind shocks. These 
shocks will encounter more cores that are now embedded in 
the larger cloud. Thus, if the cluster can trigger 14 — 25% of 
its population from the cloud cores that are active today, it 
seems plausible that it can trigger much more over the full 
lifetime of the HII region. 



with triggered stars. However, the distinct appearance of 
the two groups, and the notion of simultaneous formation 
of intermediate- mass stars and PMS siblings, favor the view 
that the star drift does not play an important role here and 
that the present locations of the two groups are directly as- 
sociated with locations of their distinct parental molecular 
cores. This is analogous to the current state of the IC 1396A 
cloud, where the two distinct stellar groups "G3" and "G4" 
inside the cloud are each associated with its own molecular 



8.3 Other aspects of triggered star formation 

8.3.1 Star Formation Efficiency 

Measured star formation efficiencies (SFEs) of entire giant 
molecular clouds ty pically vary from a few to several percent 
l|Evans et al.ll2009l . and references therein) , whereas SFEs of 
embedded clusters associated with individual massive dense 
molecular cores are significant ly higher ranging between 10% 
and 30% l|Lada fe Lada|[2003l . and references therein). 

For the total triggered population in/around IC 1396 A 
of > 140 stars, assuming an average mass of ~ 0.5 Mq, con- 
sistent with a standard IMF, and adding individual masses 
of likely triggered high-mass B-type stars in the region, 
> 10 M Q for the B3+B5 binary star HD 206081 (Chandra 
source #382) and several M for the B3 star HD 239710 
(Chandra source # 135; Figure I13t 0. we obtain an esti- 
mated total triggered stellar mass of M sta r ~ 90 Mq. The 
mass of the observed mo lecular gas in IC 1396A is only 
M oa s ~ 200 M Q today (|Patel et al l Il995l ; IWeikard et all 
Il996h . The inferred current SFE for the IC 1396 A stellar 
population is thus M s tar/(M star + M gas ) ~ 30% at the top 
of typical SFE range measured in other active star forming 
regions. 

Even higer SFEs of > 40% were recently reported in 
a number of c ases of likely trigge red star form ation, e.g., 
Cep B/OB3b (ICetman et al.ll2009ri, Sh 2-23 3IR l|Yan et all 
120101 ), L935/NGC 7000 (|Touiima et aUl201lD . However, the 
gas mass detected today could be only a fraction of the orig- 
inal cloud mass, so the SFE would be lower with respect to 
the original gas mass. Thus, the high values of apparent SFE 
in IC 1396A and the other regions of triggered star forma- 
tion indicate either an efficient star formation or an efficient 
removal of molecular gas by photoevaporation and stellar 
winds. 



8.3.2 Long-lived RDI triggering in a clumpy cloud 

Intermediate-mass stars are often born surrounded by small 
;roups of a few to a few dozen lowe r-mass PMS stars 
Testi et al.lll999l ; lAdams fc Mverdl200ll ). The stellar group 
"G2" is composed of > 15 YSOs around the B3+B5 binary 
system HD 206081 (Figure |T3t,d). The high apparent disk 
fraction of ~ 60%, the presence of 6 high-accretion YSO sys- 
tems from Bll, and the filamentary structure seen at 24/im 
coincident with the stellar group suggest that the group was 
formed very recently. The projected distance between "G2" 
and the stellar group "Gl", located right along the optical 
Rim A, is 1.5 pc. Theoretically, it is possible that stars re- 
cently formed in "Gl" having acquired velocity dispersion 
> 1 km s _1 have drifted 1 — 2 pc in all directions, popu- 
lating the field around "Gl" including the location of "G2" 



This implies continuing star formation over millions of 
years when the RDI mechanism occurs in a larger than the 
presently seen IC 1396A cloud composed of multiple clumps. 
This view is consistent with the recent global picture of an 
Hn region expanding into a non-uniform clumpy molecu- 
lar cloud, where triggered star formation would take place 
through simultaneous enhancement of density and global 
radiation-dri ven implosion of n umerous pre-existing molec- 
ular clumps l|Walch et alj|201ll ). 



9 CONCLUSIONS 

Signatures of small-scale triggered star formation on the pe- 
riphery of H n regions can be erased on short timescales 
and thus, it is often difficult to identify sites of triggered 
star formation and to measure the impact of triggering pro- 
cesses. There are only a handful of studies in the literature 
that have attempted to estimate the contribution of trig- 
gered star formation from expanding Hn regions. Infrared 
emitting protostars in cloud cores and larger populations of 
PMS stars with spatial gradients in stellar age indicative of 
long-standing triggering processes have both been seen as- 
sociated with several clouds on the periphery of expanding 
H II regions. Much of the progress is due to the fruitful com- 
bination of X-ray and optical/infrared surveys; the former 
captures the older diskless PMS population while the latter 
detects the disky and accretion population. 

IC 1396 is a nearby, large shell-like Hn region, where 
traces of recent triggered star formation are still evident. We 
present new Chandra X-ray data and auxiliary IR Spitzer 
and optical FLWO/LAICA data of YSOs associated with 
the central cluster Trumpler 37 and the adjacent bright- 
rimmed cloud IC 1396A. These data are merged with the 
data from previously published optical-IR stellar catalogs of 
the region. The goal of this work is to identify, understand, 
and quantify the triggered stellar population in the region 
associated with the IC 1396A cloud and to estimate the 
contribution of triggered star formation to the entire Hn 
region. 

Out of 415 X-ray sources detected in the 
Tr 37/IC 1396A area, 175 are identified as YSOs and 
classified on 124 diskless and 51 disky objects (*J3]). The 
majority of the remaining X-ray sources are contaminants, 
mainly extragalactic sources unrelated to the region 
In addition to the 175 X-ray emitting YSOs, we identify 
42 non-Chandra IR-excess low-mass members of the region 
O. Combining these data with the previous optical (SA05, 
Bll) and IR (MC09) catalogs of YSOs yields a total of 259 
YSOs in/around IC 1396A cloud within the ACIS field, half 
of which are new members discovered in this work (^6]). 

The age and spatial distributions of young stars reveal 
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a spatio-temporal gradient of stars from the IC 1396A cloud 
toward the primary ionizing star HD 206267 (Sj5] and [6]). 
Young stars are found to be clustered along the ionized rim, 
both outside and inside the cloud. The clustering of young 
stars in front of the cloud is consistent with the previously re- 
ported arc-shaped layer of the low-mass high accretion stars 
(Bll) and/or the ~ 1 Myr old optical stars (SA05). Quanti- 
tatively, the number of stars clustered in front of the cloud 
is consistent with the 2-fold increase in the stellar surface 
density at the western border of the central Tr 37 cluster 
(flol and |8.2.2|I . The stellar age increases from < 1 Myr in- 
side the cloud, to < 2 — 3 Myr in front of the cloud, to 
~ 4 Myr towards the central Tr 37 cluster. Combined with 
the results from the previous studies, these findings provide 
strong evidence for RDI triggered star formation over mil- 
lions of years in the Tr 37/IC 1396A region f ggTjl . Based 
on this evidence, we believe that the majority of the YSOs 
found inside the cloud and the majority of the stars found 
clustered in front of the cloud were likely formed through 
the RDI process extending over 2 — 3 Myr. 

By scaling the XLF/IMF measurements obtained for 
different YSO samples in the region (JFJ to the number of 
the observed triggered stars, we estimate a total of > 140 
triggered stars down to 0.1 Mq in/around the IC 1396A 
cloud, constituting > 25% of the total popualtion in the 
central part of the region today ( i]8.2.2|l . Furthermore, with 
a few additional assumptions we estimate the contribution 
of the triggered stellar population to the total population of 
the entire IC 1396 Hn region to be > 14% f fllTlOj) . Thus, 
the currently active clouds on the periphery of IC 1396 add 
at least 14 — 25% to the current cluster population. When 
past and future clouds are considered, the contribution of 
triggering to the star formation in the molecular cloud could 
be much higher. 

The inferred apparent value of star formation efficiency 
for the triggered population in IC 1396A is 30% using cur- 
rent gas cloud masses ( H8.3.ip . Our findings support a pic- 
ture of continuing star formation over millions of years when 
the RDI mechanism occurs in a larger than the presently 
seen IC 1396A cloud composed of multiple clumps ( H8.3.2| l. 
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